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1. INTRODUCTION AND OBJECTIVES 

1.1 Introduction 

This Final Report has been prepared by Honeywell Engines & Systems, Phoenix, Arizona, a 
unit of Honeywell International Inc., documenting work performed during the period August 1999 
through December 1999 for the National Aeronautics and Space Administration (NASA) Glenn 
Research Center, Cleveland, Ohio, under the Small Engine Technology Program, Contract No. 
NAS3-27483, Task Order 23, Wing Reflection Code. The NASA Task Monitor was Mr. Robert A. 
Golub, NASA Langley Research Center, Mail Code 461, Hampton, Virginia 23681-0001; 
telephone: (757)864-528 1 . The NASA Contract Officer was Ms. Linda M. Kendrick, NASA 
Glenn Research Center, Mail Code 500-305, Cleveland, Ohio 44135-3191; telephone: (216)433- 
2407. 

The work performed under Task 23 consisted of the development and demonstration of 
improvements for the NASA Aircraft Noise Prediction Program (ANOPP), specifically targeted to 
the modeling of engine noise enhancement due to wing reflection. 

This report focuses on development of the model and procedure to predict the effects of wing 
reflection, and the demonstration of the procedure, using a representative wing/engine 
configuration. 

1.2 Objectives 

The primary function of the ANOPP program 11, 2) * is to provide the best, currently-available 
methods to predict aircraft noise. As new methods and engine acoustic data become available, 
ANOPP prediction modules can be improved or new modules can be added. 

One aspect of noise modeling that has not previously been addressed in ANOPP is the effect 
of engine noise reflected from the wing surface, as perceived by a ground observer. Noise at the 
receiver is enhanced due to reflection of the noise incident on the wing surface. 

The objective of this task was to develop a procedure to predict the effects of wing reflection 
on farfield acoustics, and to provide a Fortran prediction code to the ANOPP noise prediction 
program for wing reflection modeling. 

1.3 Summary 

A wing reflection model was successfully developed and prepared for installation in the 
ANOPP program. The model represents the increase in noise at a ground observer position that 
results when the aft fan duct source reflects off of the wing surface. 

The degree of modeling accuracy required to predict the wing reflection behavior was 
determined using the Raynoise® [3] ray tracing program. This tool was applied to a series of 
wing/engine configurations to assess the impact of extended flaps, etc., on the noise reflected by 
the wing. The wing reflection studies performed with Raynoise® showed the importance of the 


* References are listed in Section 8. 
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wing and flap as reflecting surfaces. As a result, the wing reflection model developed for ANOPP 
included multiple reflecting surfaces, in order to adequately represent the wing geometry. 

Based on the conclusions reached in the Raynoise® study, ANOPP-compatible software was 
developed to model the wing reflection, in order to obtain more accurate farfield acoustics 
predictions. The wing reflection model was initially implemented in the Engines & Systems noise 
prediction program, GASP 141 , and was demonstrated with a typical high-wing regional transport 
aircraft. Use of the wing reflection model resulted in an increase in observed noise from the aft fan 
duct, which was consistent with the predictions from the Raynoise® analyses. 

Along with the wing reflection software, documentation was prepared, including a User Manual 
and a Technical Manual for the ANOPP wing reflection module. 
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2. TECHNICAL APPROACH 

2.1 Approach 

Experimental investigations have shown that the impact of wing reflection on farfield 
acoustics can be significant^. In particular, for regional aircraft with engines mounted under the 
wing, engine noise reflected from the wing can make a substantial contribution to the total noise 
detected by a ground observer. 

NASA’s ANOPP noise prediction program does not currently have a model to predict noise 
due to wing reflection. In order to improve the prediction capability of ANOPP, for under-the- 
wing engine mount configurations, a wing reflection model was needed. 

As a first step toward developing the wing reflection model, a study was performed to 
determine the degree of modeling accuracy required to capture the primary elements of the 
reflected noise. The Raynoise® ray-tracing program was applied to examine a series of 
wing/engine configurations, and assess the impact of extended flaps, etc. on the noise reflected by 
the wing. This effort is discussed in more detail in Section 3. 

In addition, various effects were examined using simple analytical estimates to determine 
what assumptions were justified, and whether a particular effect had to be included in the final 
reflection model. A criterion for determining the level of accuracy required involved the 
determination of the impact of the contribution from a particular source, during a 2-second 
integration of flyover noise exposure. If a particular source was not estimated to contribute 
measurably to the 2-second integrated noise level at the ground observer, then the contribution was 
considered to be secondary, and was not included in the model. Details of the model development 
are presented in Section 4. 

The general model for wing reflection that was developed, based on the conclusions reached 
in the Raynoise® studies and the additional analytical estimates, then formed the foundation of the 
wing reflection analysis procedure that was implemented in a Fortran software module compatible 
with the ANOPP program. The wing reflection module was first installed in the Engines & 
Systems noise prediction program, GASP, where a series of analyses were performed to 
demonstrate that the farfield acoustics predictions accurately modeled the wing reflection effects. 
The results were compared with Raynoise® predictions for equivalent wing/engine configurations. 
The GASP analyses are summarized in Section 5. 

After demonstrating the method in the GASP program, the wing reflection module was 
installed in ANOPP. This effort is discussed in Section 6. 

2.2 Selection of Configurations 

Relevant wing/engine configurations were required for both the Raynoise studies and the 
demonstration cases for the new wing reflection model. It was desired to use a typical regional 
configuration, with engines mounted under-the-wing. The RJ100 regional transport, shown in 
Figure 1, was selected for the wing reflection studies. It is a high-wing transport with four engines 
mounted under-the-wing. The wing has a single, unbroken, tabbed Fowler flap system, which 


3 



extends over 78 percent of the wing span. Engines & Systems’ LF507 engines are on the aircraft, 
and therefore sufficient information is available to model the aircraft and engine geometry. 

For purposes of the Raynoise' studies and new w ing reflection model predictions, the aircraft 
is modeled as simply a wing/engine combination. The fuselage and empennage are not included. 

The certification flight profiles were used for take-off, sideline, and approach flight 
conditions. Analyses were performed for configurations with and without flaps, as appropriate to 
the flight condition. The primary flight conditions and wing configurations analyzed were as 
follows: 

• Cutback Take-Off: Wing with Flaps Retracted 

• Sideline: Wing with Flaps Retracted 

• Approach: Wing with Full Flaps (33 Degrees) 



Figure 1. The RJ100 Regional Transport Served as the Basis for the Wing/Engine 
Configuration l sed in the Wing Reflection Studies. 
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3. WING REFLECTION USING RAY-TRACING PROGRAM 

3.1 Ray-Tracing Predictions 

The commercial acoustic ray-tracing program Raynoise® from LMS was used to predict the 
effects of engine noise reflections off an aircraft wing. The detailed noise contours provided by 
Raynoise® gave an understanding of both the gross effects of wing reflections as well as some of 
the nuances related to source position and wing geometry. Results from the ray-tracing program 
were used to optimize the methodology for the wing reflection model to be developed for NASA’s 
ANOPP noise prediction program. 

3.2 Description of Raynoise® Program 

Raynoise® models the physics of acoustical propagation, including specular and diffuse 
reflections against physical boundaries, wall absorption and air absorption, diffraction, and 
transmission through walls. 

Raynoise® combines a mirror image source method (MISM) with a ray tracing method 
(RTM), called the triangular beam method. Beams having a triangular cross section are emitted 
from the noise source. The number of beams the user selects determines the angular coverage of 
each beam. Having a triangular shape allows for complete coverage of the surrounding volume 
without the beams overlapping. The axis of each beam is traced while assuming spherical 
divergence and the corresponding decrease of sound intensity with distance. Absorptive properties 
of surrounding surfaces are accounted for with appropriate energy reductions with each reflection. 
Receiver locations are modeled as points and the contributions from all beams arriving at each 
receiver point are summed. Use of beams eliminates the need for time-consuming visibility testing 
inherent in the MISM while maintaining the advantage of simple receiver definition (point 
definition for MISM vs. volume required by RTM). 

Single point sources are used in Raynoise® to model compact noise sources, while an array of 
point sources is used to model a distributed source. Both incoherent and coherent sources are 
available. Only incoherent sources were used in the analyses presented here, because it was 
assumed the relative phase between direct and reflected rays is not preserved in actual atmospheric 
propagation. The user may specify the angular extent of the source propagation (i.e., spherical, 
hemispherical radiation). Raynoise® can handle both specular and diffuse reflections of multiple 
order — beams can reflect more than once. Specular reflection was assumed for all studies reported 
here. 

The geometry of a model along with the receiver locations can be constructed in Raynoise ; 
however, it was found that creation of these components using another piece of software was more 
efficient. MSC Patran® was used to manipulate IGES files for the aircraft under study and to 
generate the necessary model and receiver meshes. 

3.3 Raynoise® Models of the RJ100 

The model that was used for the Raynoise® study was derived from the RJ100 regional 
transport. The aircraft model was simplified to include only the port wing, flaps, outboard LF507 
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engine, and pylon (Figure 2). The wing and flaps were modeled as flat plates, and the engine 
nacelle and pylon surfaces were approximated by sets of planar surfaces (Figure 3). 

Two models were considered for the pressure side of the wing. Simple connection of the 
leading and trailing edges produced the first wing model, described by a single element (Figure 3). 
The second wing model was slightly more complex and used two elements to estimate the 
curvature of the wing (Figure 4). Both models were run in three different configurations: engine- 
wing, engine-wing-flaps, and engine-wing-flaps-pylon. Flaps were modeled in the fully extended 
position. 

Material properties were assigned to elements within the model. The wing and flaps were 
treated as perfectly reflective surfaces. All other surfaces were modeled as completely absorptive. 
This eliminated the possibility of reflections from the planar surfaces that represented the nacelle 
and pylon surface curvatures. 

Raynoise® was used to model fan discharge noise from the LF507 engine. Because the LF507 
installation has a separate flow nozzle, the source was placed in the fan bypass exit plane and the 
core nozzle was removed. A 4000 Hz tone was modeled to simulate the fan blade passage 

frequency. However, since an incoherent source was assumed, the predictions are valid for any 
frequency. 

Three different source configurations were used. The first was a point source located on the 
centerline of the outboard engine in the plane of the fan bypass duct exit, as is shown in Figure 3. 
The point source was configured to radiate over the aft-facing hemisphere. The second source 
model included both the inboard and outboard engines modeled using two point sources that both 
radiated only to the aft-facing hemisphere (Figure 5). The inboard point source was positioned on 
the centerline of the inboard engine in the plane of the fan bypass duct exit. A more realistic 
distributed source was also used to model the outboard engine, making up the third source model. 

Twelve point sources were aligned in a ring representing the annular fan bypass duct exit (Figure 

6 ). 


A ground plane was modeled as shown in Figure 7. The aircraft was positioned at an altitude 
of 394 feet above a plane extending 700 feet in every direction. Grid spacing on the ground was 20 
feet. 


3.4 Model of Single-Element Wing with Point Source 

As a baseline, noise contours were generated for the engine only with a single point source 
(no wing reflection), as shown in Figure 8. 

The effect of reflections off the single-element wing (no flaps or pylon) is shown in Figure 9. 
With the source located just forward of the wing leading edge (Figure 3), the boundary of the 
region affected by wing reflections was apparent in both the contour and delta plots. Noise levels 
were increased by approximately 3 dB at each receiver location in the affected region, because each 
receiver (grid point) “saw” two rays that had traveled nearly the same distance from the engine. 

Addition of the flaps increased noise levels on the ground plane, immediately aft of the 
aircraft position, as three and sometimes four rays propagated to the same location (Figure 10). 
Rays initially reflected to the aft starboard side of the ground plane were shielded by the flaps. 
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Single reflections were directed forward of the aircraft by the flaps. The widths of these forward 
regions, where noise levels increased, were 40 to 60 feet for a flyover altitude of 394 ft. 

In an attempt to determine the influence of an engine pylon on the reflection of fan exit noise 
off a wing, the geometry of the RJ 100 pylon was described using a simple four-element model 
(Figure 3). The beam method of Raynoise® breaks down for sharply curved surfaces like the 
underside of the engine pylon, thus the simple model. Because the pylon could not be accurately 
simulated, it was decided that the assignment of a completely absorptive material to the four pylon 
surfaces would be sufficient to determine the extent of the ground plane affected by this region of 
the wing. Figure 1 1 shows that, with a point source located in the lateral center of the pylon, a 
“trough of silence” was predicted behind the aircraft. This trough was approximately 160-feet 
wide. In actuality, the fan exit noise source is not a point but is rather a distributed source having a 
diameter nearly three times the width of the pylon. The results with a distributed source are given 
later. 

3.5 Model of Double-Element Wing with Point Source 

Division of the underside of the wing into two elements produced two distinct regions 
affected by reflections. Figure 12 shows that the forward portion of the wing reflected fan noise 
down to the region below and just aft of the aircraft, while reflections off the aft portion of the wing 
increased noise in the region farther behind the aircraft. A large region unaffected by reflections 
was predicted between them. Reflections were directed forward of the aircraft on the starboard 
side, due to the slope of the forward wing surface. 

Addition of flaps and an absorptive pylon produced similar trends to those with the single 
element wing. The flaps increased noise forward and just aft of the aircraft while the absorptive 
pylon created a trough that was unaffected by reflections (Figure 13). An interesting result was the 
appearance of reflections between the regions influenced by the forward and rearward wing 
elements. These were the result of double reflections off the wing and flaps and influenced isolated 
areas approximately 40-feet wide. 

3.6 Model of Inboard and Outboard Engines with Point Sources, Single Element Wing 

Inclusion of the inboard engine produced very similar noise contours to those predicted with 
only the outboard engine (Figure 14). The close proximity of the two noise sources and the wing 
relative to their altitude resulted in only slight changes in the sizes of the regions affected by wing 
reflections. A comparison was made between the noise contours predicted with both the inboard 
and outboard engines and the contours predicted with only the outboard engine having twice the 
source strength (plus 3 dB). Figure 15 shows that over most of the ground plane the two models 
gave the same results. Differences arose in narrow regions near the edges of those areas influenced 
by reflections because of the differences in source locations. 

3.7 Model of Outboard Engine with Ring Source, Single Element Wing 

The point sources used in the 12-point distributed ring source model each had the same sound 
power level as the single point source used in the models presented above. For that reason, the 
overall noise levels along the ground plane were higher, as shown in Figure 16 for the distributed 
source baseline with no wing reflection. 
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Use of a spatially distributed source served to “blur” the edges of the regions affected by wing 
reflections (Figure 17), but produced essentially the same gross behavior as seen for the point 
source. 

With flaps deployed, the area forward of the aircraft that was influenced by reflections off the 
aft flap (tab) became wider (Figure 18). Recall that the single point source was predicted to create 
a band approximately 40-feet wide (Figure 10). The ring source produced a band four times that 
wide. The radius of the distributed source was of the same order of magnitude as the chord length 
of the tab, and the two were relatively close to each other. Therefore, the range of possible 
reflection angles off the tab relative to the engine centerline was much greater than was produced 
with the single point source. A decrease was seen in the ratio of rays reflecting off the wing to 
those reflecting off the flaps and propagating to locations just aft and starboard of the aircraft. This 
is shown in Figure 18, where the maximum increase in noise due to wing-flap reflections was 
predicted to be just over 4 dB. Compare this result with the result from the single point source 
model (Figure 10), which showed a predicted increase of 6 dB at some locations. These 
predictions indicate that modeling the fan exit noise as a point source may give a conservative 
(high) estimate. 

It was seen that the presence of an engine pylon did not have a significant effect when the fan 
source was modeled as a ring of point sources with diameter approximately three times the width of 
the pylon (Figure 19). The “trough of silence” directly behind the aircraft essentially disappeared. 
The pylon provided some shielding of rays propagating from one side of the pylon to the wing on 
the opposite side, resulting in a slight depression of the noise in the region behind the aircraft. It is 
likely that in reality, reflections from the pylon would make up for the noise decreases due to 
shielding, thereby giving a result very similar to that for the model without a pylon (Figure 18). 

3.8 Effect of Double Reflections 

It was discovered that rays that had reflected twice - first off the wing and then off a flap - 
influenced noise levels over a certain area of the ground plane. The locations where this was true 
are shown in Figure 20 for the single-element wing and flaps model. A narrow strip of the ground 
plane (never more than 40-feet wide) just aft of the aircraft was affected. Extraction of the noise 
contributions of the rays that reflected twice from the total noise levels at each affected point gave 
the deltas shown in Figure 21. Double reflections were predicted to increase noise levels 1.2 dB at 
most of the locations marked in Figure 20. For an aircraft forward velocity of 300 ft/sec, the region 
influenced by double reflections would propagate over a stationary ground point in 0.13 second. 

3.9 Conclusions 

Several conclusions were reached, based on the results of the Raynoise® analyses. 

1 ) Many of the issues examined produced only localized effects, which would have minimal 
impact when integrated over a 2-second time interval: 

- Fully-extended flap surfaces 

- Engine position relative to the wing (i.e., inboard versus outboard mount locations) 

- Multiple reflections of a single ray 
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2) Use of a point source, as opposed to a distributed source, provided greater simplicity for the 
model, and produced the same general behavior for wing reflection effects, although it possibly 
caused predictions to be somewhat conservative. 

3) Modeling of the pylon was somewhat misleading, because the pylon surfaces were treated as 
completely absorptive. In reality, the reflections off the pylon surfaces would eliminate much 
of the “trough of silence”, and a realistic distributed source would also reduce the trough effect. 
Thus, eliminating the pylon from the point-source model would very likely produce more 
accurate results than attempting to model it in a simplified manner. 

4) For this aircraft configuration, use of the double-element wing showed major changes in the 
predicted noise behavior. However, modeling requirements may differ from case to case, and 
would be expected to be a function of the given wing geometry. Therefore, the accuracy of the 
noise predictions will be dependent on the establishment of best practices regarding modeling 
of the wing. Similarly, the impact of extended flaps may vary from case to case, and also 
should be examined in more detail, to establish best practices. 

5) In general, it may be observed that there is no significant benefit to analyzing a complex 
geometric model with a simple point source. Instead, application of a simple planar wing/flap 
model produces reasonable results with a point source. 

The insight into wing reflection modeling gained through the Raynoise® studies was applied 

to the development of the reflection model for the ANOPP noise prediction program. 
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Figure 2. Outboard LF507 Engine on RJ100 Aircraft. 



Figure 3. Raynoise" Model of the RJ100 Outboard Engine With Single-Element Wing. 
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Figure 4. Raynoise" Model of the RJ100 Outboard Engine With Double-Element Wing 
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Figure 5. Raynoise" Model of the RJ100 Inboard and Outboard Engines, With Single 
Element Wing. 
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Figure 6. Raynoise' Model of the RJ100 Outboard Engine Using a Distributed Source, W ith 
Single-Element Wing. 



Figure 7. Ground Plane for All Predictions, Aircraft at Altitude of 394 Ft. 
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Figure 8. Noise C ontours for Engine-Only Configuration, With Point Source. 
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Figure 9. 


Engine With Single-Element W ing: Noise Contours, Noise Deltas From Engine- 
Only Configuration. 
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Figure 10. Engine With Single-Element Wing and Flaps: Noise Contours, Noise Deltas From 
Engine-Only Configuration. 
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Figure 1 1. Engine With Single-Element Wing, Flaps, and Absorptive Pylon: Noise Contours, 
Noise Deltas from Engine-Only Configuration. 
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Figure 12. Engine With Double-Element W ing: Noise Contours, Noise Deltas From Engine- 
Only Configuration. 
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Figure 13. Engine With Double-Element Wing, Flaps, and Absorptive Pylon: Noise 
Contours, Noise Deltas From Engine-Only Configuration. 
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Figure 14. Inboard and Outboard Engines With Single-Element Wing, Flaps: Noise 
Contours, Noise Deltas From Engine-Only Configuration. 
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Figure 16. Noise Contours for Engine Only, 12-Point Ring Source. 
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Figure 17. Outboard Engine, 12-Point Source, Single-Element Wing: Noise Contours, Noise 
Deltas from Engine-Only Configuration. 
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Figure 18. Outboard Engine, 12-Point Source, Single-Element Wing and Flaps: Noise 
Contours, Noise Deltas From Engine-Only Configuration. 
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Figure 19. Outboard Engine, 12-Point Source, Single-Element Wing, Flaps, Absorptive 
Pylon: Noise Contours, Noise Deltas From Engine-Only Configuration. 
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Figure 20. Double Reflections and the Locations Influenced by Them, Single Element Wing 
and Flaps, Point Source. 
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Figure 21. Comparison of Wing-Flap Reflections With and Without Double Reflections, 
Single-Element Wing and Flaps, Point Source. 



4. DESCRIPTION OF WING REFLECTION MODEL FOR ANOPP 

4.1 Aspects of Noise Modeling for Engine Under-the-Wing Configurations 

As shown in Figure 22, wing reflection was but one of many issues which needed to be 
addressed, in developing the prediction module for an engine under-the-wing configuration, 
including the noise source, the directly-radiated wave, the reflected wave, entrainment-induced 
noise, and various effects on noise propagation. Each of these will be discussed in more detail in 
the following sections. 

4.1.1 Noise Source 

The noise source was limited to fan bypass duct noise, which is the main contributor to wing 
reflection noise. A point source was assumed, with a spherical sound wave radiating in an aft- 
facing hemisphere, as shown in Figure 23. This was consistent with Heidmann's 161 observation that 
the fan discharge duct broadband noise at 90 degrees was reduced approximately 5 dJB from its 
maximum level, and the contribution at angles less than 90 degrees (forward of the duct exit) was 
negligible. 

The sound field was assumed to be axisymmetric, in the azimuthal direction; however, polar 
directivity was accounted for in the model. 

4.1.2 Directly-Radiated Wave 

In the wing reflection model, the directly-radiated wave was treated in a manner consistent 
with the current ANOPP program. All of the significant effects were accounted for, including 
directivity, atmospheric attenuation, spherical spreading between the source location and the 
ground observer location, ground reflection, and extra-ground attenuation (for sideline predictions). 

4.1.3 Reflected Wave 

Among the modeling considerations which had to be taken into account for the reflected wave 
were the nature of the reflecting surfaces, the directivity of the source, and the issue of single 
versus multiple reflections of a single ray. 

The aircraft configuration itself had to be modeled, as part of the definition of reflecting 
surfaces. For the ANOPP model, it was assumed that reflection could occur only from the 
wing/flap system. There would be no contribution from the fuselage, empennage, or engine pylon 
surfaces. Influence of the pylon was examined during the Raynoise® studies, and it was concluded 
that the pylon could be neglected in the final wing reflection model. 

A maximum of three wing/flap panels were included in the model. The flap system was 
assumed to consist of continuous spanwise flap panels on each wing, with no jet efflux cut-outs. 
Tabbed Fowler flaps could be modeled, including gaps along the hinge centerlines. 

The effect of leading edge slats on the wing was neglected. The slats would have a secondary 
effect as reflecting surfaces, relative to the wing/flap system. In addition, the assumption of an aft- 
facing hemispherical source model would limit the range of directivity angles at which a slat 
reflection would occur. The size of the slats would also limit the time interval during which 
flyover noise reflected from the slat would reach the receiver. 
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The characteristics of the reflecting surfaces also had to be considered, when developing the 
wing reflection model. The wing and flap panels were modeled as planar reflecting surfaces. Each 
reflecting panel was assumed to be an infinitely hard acoustic surface (i.e., perfectly reflecting). 

Directivity of the source noise had to be included in the wing reflection model, because the 
polar angle at which the reflected ray leaves the engine may vary significantly from that of the 
direct ray, as shown in Figure 24. The difference in angle would depend on the position of the 
reflecting surface relative to the source, and the position of the source relative to the receiver. This 
difference in angle would result in different source strengths for the direct and reflected rays. 

The issue of single versus multiple reflections was considered during the Raynoise® studies. 

It was determined that the impact of multiple reflections of a single ray was negligible, for the case 
under study. As a result, the ANOPP wing reflection model assumed that a single ray would reflect 
from at most one surface (wing, flap, or flap tab), and reflections of a single ray off of multiple 
surfaces would be neglected. In addition to the Raynoise studies, a similar case (RJ100 wing with 
LF507 engine) was modeled for a two-dimensional wing/flap system with fully-deployed flaps, as 
shown in Figure 25. A geometrical study of the reflected rays showed that instances of multiple 
reflections would occur for a range of approximately 2 degrees, in terms of aircraft-observer angle. 
This is small compared to the changes in angle during each Vi second interval of flyover noise 
predictions. Therefore, neglecting the multiple reflections should have no significant impact on the 
noise predictions. 

4.1.4 Entrainment-Induced Noise 

Entrainment-induced noise represents an additional noise source for engine under-the-wing 
configurations. This phenomenon occurs when flow is entrained between the nacelle and the wing, 
due to jet exhaust velocity. The entrainment-induced boundary layer on the wing surface is 
responsible for the noise. Entrainment effects are a function of the jet exhaust velocity, the axial 
wing position relative to the engine exhaust, and the normal distance between the wing and engine 
centerline. 

According to Wang [5] , entrainment-induced noise predominates at low frequencies. In scaled- 
nozzle tests, Wang found that the frequencies were below 500 Hz. When the data are scaled up for 
typical aircraft configurations, the frequencies at which entrainment has an impact are lowered 
(below approximately 170 Hz, for the RJ 100 with LF507 engines). Therefore, the impact is at low 
enough frequencies to be of no consequence at full scale. 

In addition, this noise source is distinct from the direct and reflected-wave noise, and 
therefore is beyond the scope of the current modeling effort. 

4.1.5 Effects on Propagation 

Propagation of the reflected wave is influenced by the same effects as the directly radiated 
wave, i.e., spherical spreading, atmospheric attenuation, and ground reflection. In the ANOPP 
model, these effects were treated for the reflected wave in the same manner as they were for the 
direct wave. 
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In addition, other effects had to be considered, which are more specific to the reflected wave. 
These included noise transmission through the jet, one aspect of which is refraction of the reflected 
wave through the jet, and interference between the directly radiated and reflected waves. 

Some of the noise reflected from the wing is transmitted through the jet as it propagates to the 
receiver. During transmission through the jet, various frequency components of reflected noise are 
either absorbed, reflected, or refracted by the jet. As noted by Wang, high-frequency components 
of noise (wavelengths small compared to the dimension of the jet) are most susceptible. Therefore, 
high-frequency sound waves would experience more attenuation when transmitted through the jet 
than when their paths do not pass through the jet. However, these effects are outside the scope of 
the current wing reflection model, and will not be addressed. 

One aspect of noise transmission through the jet, that of refraction, can be assessed in an 
approximate fashion to determine whether the impact would be significant. A reflected wave 
passing through the jet undergoes refraction, which alters the directivity and distance that the 
reflected wave travels (Figure 26). The jet exhaust velocity changes the direction of the ray, and 
the jet temperature alters the local speed of sound. For a typical flyover case, this change in 
distance compared to the overall distance traversed by the reflected wave was negligible. In 
addition, the shift in ray position was negligible compared to the distance traveled by the source in 
a Vi second interval of the flyover. Therefore, it may be assumed that refraction effects are 
negligible, in terms of distance and directivity, in the wing reflection model. 

Interference can occur due to different path lengths between the directly radiated and reflected 
noise. However, the coherence coefficient for a typical case with the RJ100/LF507 configuration 
was estimated to be 0.0257, which implies that the uncorrelated energy is within approximately 
0.1 1 dB of the initial acoustic energy. Therefore, it may be assumed that interference effects are 
essentially negligible. In addition, no evidence of attenuation due to interference was seen in the 
jet noise test data of Wang. Thus, the effect was not modeled in ANOPP. 

4.2 Implementation in Wing Reflection Algorithm 

4.2.1 Procedure for Calculation of Wing Reflection 

Based on the modeling assumptions discussed above, the wing reflection algorithm for the 
ANOPP noise prediction program was developed. The primary steps in the procedure are 
presented in the flowchart of Figure 27, and include the following: 

. Define the geometry for each reflecting panel (wing, flap, flap tab) 

. Create the image source for each reflecting panel, at each flyover location 
• Determine if a reflection will occur for a given panel 
. Compute the path length of the reflected ray, for a given panel 
. Obtain the source strength for the reflected ray, based on directivity 
. Compute the strength at the receiver, by applying appropriate corrections 
. Combine the contributions from the reflected ray(s) and the direct ray 

Each of these steps is discussed in more detail in the following sections. 
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4.2.2 Geometry Definition for Reflecting Panels 

The reflecting panels that represent the wing/flap system are modeled as planar surfaces. The 
geometry definition for each panel consists of the specification of the Cartesian coordinates (X, Y, 
Z) of each comer point of the panel; i.e., coordinates are defined for the Root Leading Edge 
(RLE), Root Trailing Edge (RTE), Tip Leading Edge (TLE), and Tip Trailing Edge (TTE), as 
shown in Figure 28. The coordinates of the Tip Trailing Edge are used only to define the wing 
boundary. The plane of the wing is defined by the three points: RLE, RTE, and TLE. Any twist in 
the wing is ignored. 

The Cartesian coordinates of the wing boundary points are defined relative to the local origin 
of the engine centerline at the fan bypass duct exit plane. Local axes are shown in Figure 28. 

Then, the engine fan duct exit and wing coordinates are transformed into a global coordinate 
system consistent with the receiver location on the ground (Point R). This transformation must 
take into account the aircraft attitude and position at the particular time of the observation. 
Therefore, the global coordinates must be calculated at every flyover position. 

Three wing/flap panels are currently permitted in the model. This allows the description of a 
main wing panel, a flap panel, and a flap tab panel. However, if flaps are not included in the 
model, or are considered to have negligible impact on reflection compared to the wing itself, the 
wing may be modeled using multiple panels, instead. 

4.2.3 Creation of Image Source 

The reflection model uses an image source point, which is positioned on a normal to the 
reflecting plane, at a distance equivalent to that between the source and the reflecting plane, but 
located on the opposite side of the reflecting plane. The orientation of the image source is shown 
in Figure 29. 

To obtain the image source point, a normal to the reflecting plane is created, which intersects 
the original source point (S). Then Point IS, the image source, is located on the normal line. Points 
S and IS are equidistant from the reflecting panel. Note that the normal line between Points S and 
IS does not have to intersect the wing panel itself, only the plane extending through the wing panel. 

First, obtain the components of a vector normal to the reflecting surface, by computing the 
cross-product of vectors along the wing root and leading edge: 


where: 


Ax = Av * A 7 rr — Ay. r * A 7 

"root J LE ^*~root 

A >’„ = *^LE~ A * LE * ^ ™ , ) 

Az„ = A* ™ , * Ay^ - Ax^ * Ay rool 

A * root ~ X RTE ~ X RLE 
by root “ y RTE ~ y RLE 
^Zroot ~ Z RTE ~ Z&LE 
A * LE = X TLE ~ X RUE 
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^LE ~ yUE y RLE 

AZle = ^tle ~Zrle 

Next, the coordinates of a point, N, on a line normal to the reflecting surface and passing 
through the engine source point, S, are obtained as follows: 


X N ~ X S + 

y N =y s +A ^ 

z N — z s + A z N 

Then, the coordinates of the intersection of the line S-N with the plane of the reflecting 
surface are computed. The intersection point is identified as Point L The coordinates of the 
intersection point are determined by solving a set of three equations in three unknowns (*/> y,> 
and z , )- Two of the equations are produced by the two-point form of the equation for the line S-N: 


y/-ys _ 0 

x n~ x s y N -y s 



The other equation comes from the three-point form of the equation for the reflecting plane: 


X I ~ X RLE 
X RTE X RLE 
X TLE ~ X RLE 


>v ->w 

Zl ~ ZrLE 

y RTE ~ y RLE 

ZRTE ~ Z RLE 

y TLE ~~ y RLE 

Ztle ~ Zrle 


Once the coordinates of the intersection point have been determined, the components of the 
vector between the source point, S, and the intersection point, I, may be computed: 


Ax, = x, -x s 

Ay, =y,-y s 
A z,=z,~ z s 

Finally, the image source point, IS, may be positioned on the normal to the reflecting surface, 
through the source point, S, by adding the above vector components to the coordinates of the Point 

I: 


x is = x i +Ax, 

y/5 = y/+Ay, 
z,s =z, +Az, 
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4.2.4 Will Reflection Occur? 

Once the image source is located, it is necessary to determine if a reflection will actually 
occur for the reflecting panel under consideration. A reflection will occur if the line passing 
between the receiver and the image source at a given flyover location intersects the plane of the 
reflecting panel within the boundaries of the actual panel, as shown in Figure 30. The orientation 
of image source and receiver changes at each interval during the flyover, and therefore, must be 
recomputed for each flyover position. 

To determine if a reflection will occur, the coordinates of the reflection point (W) on the plane 
of the wing panel must first be obtained, by determining the intersection point of the line formed by 
IS-R and the plane of the wing panel. The coordinates of the intersection point, W, are determined 
by solving a set of three equations in three unknowns ( x w , y w , and z w ). Two of the equations are 
produced by the two-point form of the equation for the line IS-R: 

*w- x r yw~y R _ Q 

x ,s~ x r y, s -yR 

X W X R ~ £ R _ q 

X I5 ~ X R Z 1 S ~ Zr 

The other equation comes from the three-point form of the equation for the reflecting plane: 


X W ~ X RLE 

y w ■ 

y RLE 

~ Zrle 


X FJE ~~ X RLE 

y RTE 

y RLE 

Zrte — Zrle 

= 0 

X TLE X RLE 

y-TLE 

— y RLE 

Z-TLE —Z-RLE 



After the coordinates of the reflection point, W, have been determined, the point must be 
tested to establish whether or not it lies within the boundaries of the reflecting plane. To 
accomplish this, the point on each wing boundary that is nearest to Point W must be located. Each 
of these points (Points Wle, W tc , and Wnp) is computed by solving a set of three equations in three 
unknowns (e.g., x w , yw rr » z w TJ r f° r the leading edge boundary). The equations are 

obtained by imposing the following conditions, e.g. for the leading edge boundary: 

1) The line W-Wle must be perpendicular to the wing boundary. This condition is represented by 
setting the dot product of the line W-Wle vector and the wing boundary line vector equal to 
zero, e.g.: 

{x w — X w ^ Rig — X TLE )+ (y w — y Wl£ \y Ri£ ~ y-TLE ) + kw ~ Z Wle \zrjx ~ Zjiz ) = 0 

2) The point Wle must lie on the wing boundary. This condition is met when the coordinates of 
the point Wle satisfy the two-point equation of the line representing the wing boundary edge, 
e.g.: 

~ x rle yy/it - y_RLE 

X TLE ~ x RLE yTLE~yRLE 
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X W LE ~ X RLE ZWiE ~ Z RLE 
X TLE “ X £L£ -ri£ “ Z RLE 

The coordinates of the reflection point, W, are then tested against the coordinates of the three 
nearest boundary points to determine if Point W is actually located on the wing or flap reflecting 
surface. 

If the reflection point W is contained within the wing panel boundaries, the reflection 
computation may continue. If Point W is outside the wing panel boundaries, no contribution to 
reflected noise occurs, and the algorithm proceeds to consider the next reflecting panel. 

4.2.5 Computation of Reflected Ray Path Length 

If a reflection has occurred for a particular wing panel, the path length of the reflected ray 
must be determined. The path length is the combined distance from the source to the reflection 
point, S-W, and from the reflection point to the receiver, W-R, as shown in Figure 3 1 . This is 
equivalent to the distance from the image source point, IS, to the receiver point, R. 

Thus, the reflected ray path length may be computed from: 


d,=V( AOMAy,) 2 +(a0 2 

where: 

A* r = X !S ~ X R 

Ay r = y, s - y R 

te r =z ls -z R 

4.2.6 Determination of Source Strength for Reflected Ray 

The strength of the point source, S, is a function of the directivity angle of the ray radiating 
from it. As illustrated in Figure 32, the directivity angles of the direct ray and the reflected ray 
could be different. Therefore, the source strength of the reflected ray must be computed separately 
from that of the direct ray. 

To obtain the source strength, the directivity angle of the reflected ray must be computed in 
three-dimensional space. The directivity angle is equivalent to the 3-D included angle between the 
engine exhaust centerline and the reflected ray. Then, the reflected ray source strength is 
interpolated from an array of source strength versus directivity angle. 

4.2.7 Computation of Strength at Receiver 

To obtain the strength at the ground receiver, the source strength of the reflected ray must be 
corrected for the same propagation effects as the directly radiated ray, i.e., atmospheric attenuation, 
spherical spreading, and ground reflection. These are all functions of the path length of the 
reflected ray, and thus will have different values than for the direct ray. Existing algorithms are 
available in the ANOPP program to compute these correction values. 
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4.2.8 Combination of Contributions from Direct and Reflected Rays 

After the corrected strength of each of the reflected rays has been computed, they must be 
combined with the strength of the direct ray, in order to obtain the total Sound Pressure Level 
(SPL) at the receiver. 

The combined contribution of rays from each of the reflecting panels is computed using the 
following equation, for every 1 /3-octave frequency band, k: 

\PA\EL 

Olog I'O 

t PANELS 

The total contribution from reflected rays may then be combined with the direct ray 
contribution, to yield the total SPL at the receiver. For every 1/3-octave band, the total SPL is 
given by: 

SPL k TOT = 10 log[lO (s/>i '-'"" l0) + io (iy>i ‘ Tor l0) ] 

4.3 Conclusion 

The wing reflection model to be applied in the ANOPP noise prediction program was 
formulated based on the results of the Raynoise ® studies and additional analytical estimates. 
Reasonable assumptions were applied, and their impact on modeling accuracy was assessed. 

The wing reflection algorithm described in this section was designed to be compatible with 
NASA’s ANOPP noise prediction program. However, rather than installing it directly in ANOPP, 
a more expedient process was to first install the wing reflection module in Engines & Systems’ 
noise prediction program, GASP. This permitted more efficient check-out and demonstration of 
the module. Then, after verifying its capability, the module was installed in ANOPP. 


Entrainment-Induced Noise 


Reflected Wave 



Figure 22. Issues To Be Considered for Noise Prediction With Engine Mounted Under-the- 
Wing. 
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Figure 23. Point Source Model for Fan Bypass Noise at Fan Duct Exit. 



Figure 24. The Directivity Angle of the Reflected Ray May Differ From That of the Directly- 
Radiated Ray. 



Figure 25. A 2-D Model of a Typical Wing/Flap System, With Fully-Extended Flaps, 
Illustrates the Occurrence of Multiple Reflections From a Single Ray. 



Figure 26. A Reflected Ray Experiences Refraction As It Passes Through the Exhaust Jet. 
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Define Geometry for Each Reflecting Panel 



Compute Reflection Point, W 




Compute Path Length of Reflected Ray 


Compute Directivity of Reflected Ray 


Obtain Source Strength of Reflected Ray 


Compute Strength at Receiver by Applying 
Appropriate Corrections 


Figure 27. Flowchart for the Wing Reflection Algorithm. 
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Figure 28. The Reflecting Panel Geometry Is Defined in Terms of Panel Boundary Points. 


IS 



I 

s 


Figure 29. Orientation of the Image Source Point Relative to the Original Source and the 
Reflecting Panel. 
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Plane of Reflecting Surface 



Figure 30. To Determine if a Reflection Will Occur, the Location of the Reflection Point Y\ 
Must Be Determined Relative to the Wing Panel. 



Figure 31. The Path Length of the Reflected Ray Is Equivalent to the Distance Between the 
Receiver, R, and the Image Source, IS. 
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Figure 32. The Source Strength of the Reflected Ray Is a Function of the Directivity Angle, 
Which Can Differ From that of the Direct Ray. 
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5. DEMONSTRATION OF THE NEW MODEL IN GASP 

5.1 Configurations and Operating Conditions 

In order to demonstrate the new wing reflection model in GASP, a series of test cases was 
performed using the wing/engine configuration of the RJ100 regional transport with LF507 
engines, as described in paragraph 2.2. For purposes of the GASP analyses, only the geometry of 
the wing/flap system and the position of the noise source were modeled. The wing and flap panels 
were defined using the four comer points on each panel, and the noise source was located on the 
outboard engine centerline at the aft fan duct exit plane. 

The GASP analyses were performed at the operating conditions for the three FAA 
certification points: 

. Cutback Takeoff 
. Sideline 
• Approach 

The analyses were performed in three modes for each of the three operating conditions: 

. No wing reflection 

. Wing reflection with the wing panel only (i.e., flaps retracted) 

. Wing reflection with the wing, flap, and flap tab panels (full flaps - 33 degrees) 

Although full flaps are only appropriate for approach, and not for takeoff and sideline 
conditions, the full flap cases were analyzed for all points, in order to obtain an estimate of 
sensitivity to flap extension. 

5.2 Results 

For each operating condition, results of the three analyses were compared by computing deltas 
of EPNL and maximum PNLTC between the wing reflection model and no wing reflection model. 
Deltas were considered for both aft fan duct noise alone and total overall noise. 

At the cutback takeoff condition, shown in Figure 33, the EPNL contribution from the aft fan 
duct increases by 2.4 dB when reflection from the wing panel (1 panel) is included. Full flap 
extension (three panels) has little impact, increasing the aft fan duct EPNL by only an additional 
0.1 dB. Overall EPNL is only increased by 0.4 dB with the wing reflection option, and there is no 
additional impact from the flaps and flap tabs. 

Figure 34 presents the results at the sideline condition. The aft fan duct contribution to EPNL 
increases by 1 .9 dB when the wing panel reflection is included. Addition of the extended flap 
panels produces another 0.2 dB EPNL. For sideline, the wing reflection effect on total EPNL is 
minimal, with a prediction of 0. 1 dB increase for wing reflection, and an additional 0.1 dB increase 
with full flaps. 

At the approach condition, the EPNL increase in aft fan duct noise due to reflection from the 
single wing panel (shown in Figure 35) is 1.9 dB, and no additional increase is seen when the fully— 
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extended flap panels are included. Overall EPNL increases only 0.5 dB with wing reflection, and 
no additional increase is seen with flap extension. Examination of individual points along the 
approach path reveals local differences in wing-only and extended-flap reflected noise 
contributions; however, the integrated effect is negligible. 

As would be expected, the maximum PNLTC deltas for each of the operating conditions are 
larger than the EPNL deltas, as shown in Figures 33 through 35. These maximum PNLTC values 
for aft fan duct noise range from 2.9 dB at approach to 3.5 dB at sideline, for reflection from the 
wing panel only. A maximum additional contribution of 0.2 dB with fully extended flaps is seen at 
cutback takeoff. The overall maximum PNLTC deltas range from 0.6 dB at sideline to 1 .6 dB at 
approach. 

It may be observed that the impact on maximum PNLTC is reasonable, representing an 
approximate doubling of the noise, due to the contribution of a reflected wave from the wing 
surface. Also, it may be noted that the impact of the flap and flap tab surfaces is not substantial, 
producing at most a 0. 1 dB increase in total noise contributed from the wing reflection model. This 
appears quite reasonable, given the size of the flap reflecting surfaces relative to that of the wing 
surface. This result is also supported by the Raynoise® analyses, which showed only a small, 
localized effect from the flaps. 

The effect of wing reflection on PNL for the approach flight profile is shown in Figure 36. In 
this case, the wing/flap system is modeled with fully extended flaps. The wing reflection 
contribution is evident after the aircraft has passed over the observer’s position. The maximum 
increase in aft fan duct PNL occurs at an engine-observer angle of approximately 1 10 degrees, and 
at approximately 125 degrees for total PNL. 

5.3 Conclusions 

The results from the new wing reflection model are generally consistent with those of the 
Raynoise* analyses. However, these wing reflection results are a strong function of the 
wing/engine configuration analyzed. As was seen in the Raynoise* analyses (Section 3), different 
relative placement of the noise source and the reflecting surface boundaries can result in different 
contributions from wing reflection. Therefore, to more fully exercise the wing reflection module, 
and to establish best practices for its use, further analysis with other wing/engine configurations is 
recommended. 
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Figure 33. Effect of Wing Reflection at Cutback Takeoff Conditions for the Wing/Engine 
Configuration of the Regional Transport Test Case. 
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Figure 34. Effect of Wing Reflection at Sideline Conditions for the Wing/Engine 
Configuration of the Regional Transport Test Case. 
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Figure 35. Effect of Wing Reflection at Approach Conditions for the Wing/Engine 
Configuration of the Regional Transport Test Case. 
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Figure 36. Impact of Wing Reflection on Noise for the Approach Flight Profile With the 
Wing/Engine Configuration of the Regional Transport Test Case. 
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6. APPLICATION IN THE ANOPP PROGRAM 


After demonstration of the wing reflection model in the GASP program, the algorithm was 
then prepared for installation in the ANOPP program as the Wing Geometric Effects module, 
which is accessed with the “EXECUTE WING” command. 

The Wing Geometric Effects module was designed to perform calculations using just the 
output from the Geometry (GEO) and Propagation (PRO) modules. The method for wing 
reflection required computing a directivity angle and adjusting the source noise level for the 
reflected ray. To implement this calculation, the noise level for the reflected ray was computed by 
interpolating the PRO received noise data table at the appropriate directivity angle. Then, a 
spherical spreading correction was applied to adjust the level for the appropriate path length. It 
was assumed that the original ground effects and atmospheric absorption calculations for the 
interpolated direct ray (at the reflected ray directivity angle) were valid for the reflected ray. 

Following installation of the wing reflection module in ANOPP, an approach analysis was 
performed, using the RJIOO regional transport, to obtain predicted farfield noise levels with and 
without the wing reflection option. A plot of flyover noise levels comparing the two analyses is 
shown in Figure 37. These results were obtained directly from the Wing Geometric Effects module 
test case. 

The ANOPP Theoretical Manual for the Wing Reflection module is included in Appendix I of 
this report. The ANOPP Wing Reflection User’s Manual is contained in Appendix II, and the 
Wing Reflection Test Case Inputs and Outputs are in Appendix III. 



Figure 37. ANOPP Predictions of Flyover Noise Levels With and Without the Wing 
Reflection Option. 
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7. CONCLUSIONS AND RECOMMENDATIONS 


Engines & Systems has completed Task Order 23, Wing Reflection Code, under the National 
Aeronautics and Space Administration (NASA)-sponsored Small Engine Technology (SET) 
Program, Contract No. NAS3-27483. This task focused on improving the engine noise prediction 
capabilities of the NASA ANOPP program through the introduction of a wing reflection model. 

The wing reflection model accounts for flyover noise enhancement due to reflection of the aft 
fan duct noise off of the wing/flap system, for under-the-wing engine mount configurations. 

The ray-tracing program Raynoise® was used to examine wing reflection behavior. Based on 
input from the Raynoise study, the wing reflection model was successfully developed and 
installed in Engines & Systems' GASP noise prediction program. The model was then 
demonstrated for a typical regional transport with engines mounted under-the-wing. As expected, 
use of the wing reflection model resulted in enhancement of the aft fan duct noise, and thus the 
overall aircraft noise, at a ground receiver location, due to reflection from the wing/flap system. 

The software module containing the wing reflection algorithm was then installed in the 
ANOPP program. In addition, ANOPP documentation for the wing reflection module was 
generated, including a Technical Manual, a User’s Manual, and a listing of input and output for a 
test case. 

Further analyses using the wing reflection module in ANOPP would be recommended, with 
different wing/engine configurations, in order to more fully exercise the model and determine its 
sensitivity to the relative placement of noise sources and reflecting surface boundaries. Also, 
further examination of the impact of flap system surfaces would be of value, to establish best 
practices for modeling small reflecting surfaces. 
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APPENDIX I 
ANOPP 

WING GEOMETRIC EFFECTS MODULE 
THEORETICAL MANUAL 


I-i 



INTRODUCTION 


The Wing Geometric Effects Module computes the effects of wing shielding and reflection 
on the propagation of noise from the engine. The wing shielding model employs the Fresnel 
diffraction theory for a semi-infinite barrier, as described in Beranek (1) and Maekawa (2), with 
modifications to treat the finite barrier presented by the aircraft wing. In the wing reflection model, 
the directly radiated wave is treated in a manner consistent with the current ANOPP program. For 
the reflected wave, the nature of the reflecting surfaces and the directivity of the source are 
modeled as well as atmospheric absorption, ground effects, and source-to-observer geometry. It is 
assumed that reflection can occur only from the wing/flap system. A maximum of three wing/flap 
panels is included in the model. The flap system consists of continuous spanwise flap panels on 
each wing, with no jet efflux cut-outs. Tabbed Fowler flaps can be modeled, including gaps along 
the hinge centerlines. The effect of leading edge slats on the wing is neglected. The wing and flap 
panels are modeled as planar perfectly reflecting surfaces. 



SYMBOLS 

A 

attenuation, dB 

Coo 

ambient speed of sound, m/s (ft/s) 

d 

distance, m (ft) 

f 

frequency, Hz 

N 

Fresnel number 

<p 2 >* 

mean-square acoustic pressure, re p 2 c* 

x, y, z 

coordinate locations 


GREEK 

A 

difference in source-receiver path length 

PoO 

ambient density, kg/m 3 (slug/ft 3 ) 


SUPERSCRIPT 

* 

dimensionless quantity 


SUBSCRIPT 

F 

flap 

1 

source location 

0 

observer location 

r 

reference standard sea level 

RLE 

root leading edge 

RTE 

root trailing edge 

T 

flap tab 

TLE 

tip leading edge 

TOT 

total 

TTE 

tip trailing edge 

W 

wing 
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INPUT 


The values of the wing (and the flap and flap tab for reflection) coordinates are provided by 
the user. The source-to-observer geometry is provided by the Geometry (GEO) Module and the 
one-third octave band noise levels being propagated to the observer are provided by the 
Propagation (PRO) Module. The frequency array establishes the independent variable values for 
the output table. 


(XRLE, yRLE, Zrle) 
(XRTEi yRTE, ZRTe) 
(XtLE, yTLE, Zjle) 
(XlTE, y-TTE, Ztte) 


coordinates of root leading edge 
coordinates of root trailing edge 
coordinates of tip leading edge 
coordinates of tip trailing edge 


Independent Variable Array 
f frequency, Hz 

Received Noise Data Table 

f frequency, Hz 

t reception time, s 

o observer index 

4c 

c a (o) speed of sound at the observer, re c r 

4c 

p a (o) air density at the observer, re p r 
<p (f,t,o)>* mean square acoustic pressure, re p^c m 

OUTPUT 

The output of this module is a table of the mean-square acoustic pressure as a function of 
frequency, reception time, and observer index corrected for wing geometry effects. 

Attenuated (or Enhanced) Received Noise Data Table 

f frequency, Hz 

t reception time, s 

o observer index 

c* (o) speed of sound at the observer, re c r 

p a (o) air density at the observer, re p r 
<p 2 (f,t,o)>* mean square acoustic pressure, re p 2 w c 4 x 
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METHOD 


Wing Geometry 

The wing (flap or flap tab) configuration is described in a local coordinate system with the 
origin positioned at the engine inlet (Point 1), as shown in Figure 1 for wing shielding. The origin 
of the local coordinate system for wing reflection is at the engine exhaust as shown in Figure 
The local coordinate system origin is assumed to be the aircraft location specified in the body 
coordinate system (see the Geometry Module) for the propagation calculation. The user must 
specify the coordinates at the wing root leading edge, root trailing edge, tip leading edge, and tip 
trailing edge, relative to the location ot the local coordinate system origin. 

Then, the local origin and wing coordinates are transformed into a global coordinate system 
consistent with the observer location on the ground (Point O). This transformation must take into 
account the aircraft attitude and position at the particular time of the observation. 



Figure 1. Coordinate System That Defines the Wing Geometry and Typical Sound 
Propagation Vector For Wing Shielding. 
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RLE 



Figure 2. The Reflecting Panel Geometry Is Defined in Terms of Panel Boundary Points 
Similar to Wing Shielding. 

Wing Shielding Model 

The location of the point representing the intersection of the line between the engine inlet 
(Point 1 ) and the observer on the ground (Point O) with the plane of the wing must be computed. 
Figure 3 illustrates the configuration of line 1-0 and the wing plane, with the intersection point 
(Point l ). The coordinates of the intersection point are determined by solving a set of three 
equations in three unknowns ( x { , >7 , and zj ). Two equations are produced by the 2-point form 
for the equation for the line 1 - 0 : 


£/ - £0 _ y_i - yp 

*1 -X 0 V’l -Vo 


( 1 ) 


x i~ x o -1 ~ z o 

The other equation comes from the 3-point form of the equation for the wing plane: 

x ( ~ x RLE y l RLE - / - - RLE 

x RTE ~ x RLE y RTE ~ RLE Z RTE ~ - RLE = 0 
X TLE ~ x RLE >' TLE ~ >' RLE - TLE ~ - RLE 
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Because four points have been specified to describe the boundaries of the wing, the wing 
surface may not actually be planar. However, for the purpose of determining the intersection Point 
/ , the assumption is made that the wing plane is described by the points at the root leading and 
trailing edges, and the tip leading edge. The intersection point (Point I ) may or may not be 
located within the boundaries of the wing surface. 

Now, the point on each wing boundary which is nearest to Point I must be located, as shown 
in Figure 3. Each of these points (Points W^, Wte, and Wnp) is computed by solving a set of 
three equations in three unknowns (e.g., x W[£ , yw LE > and )• The equations are obtained by 

imposing the following conditions: 

1) The line I -W must be perpendicular to the wing boundary. This condition is represented 
by setting the dot product of the line I -W vector and the wing boundary line vector equal 
to zero, e.g.: 

(*/ - X W[£ )f, x RLE - X TLE ) + (?/“ ywur \yRLE~ yTLE ) + ( Z I ~ ){ Z RLE ~ Z TLE ) = 0 

2) The point W must lie on the wing boundary. This condition is met when the coordinates 
of the point W satisfy the two-point equation of the line representing the wing boundary 
edge, e.g.: 

X W LE ~ X RLE yw LE ~ yRLE _ Q (5) 

X TLE ~ X RLE yTLE ~ yRLE 

X WLE ~ ^le z W le ~ ZRLE _ Q (6) 

X TLE ~ X RLE Z TLE ~ Z RLE 

It is necessary then to determine if the intersection point / actually is located within the 
boundaries of the wing. If it is outside the wing, then no attenuation of the noise source is present. 
However, if Point / lies on the wing surface, then the Fresnel diffraction theory is applied to 
determine the level of attenuation. 

Assuming that Point / is located within the boundaries of the wing, then the attenuation of 
the noise source due to wing shielding must be determined for each diffraction edge (i.e., wing 
boundary edge). For each diffraction edge, three distances must be computed, as shown in Figure 
4: 


1) The direct source-receiver path length, from Point 1 to Point O, dio, 

2) The distance from Point 1 to the closest point on the diffraction edge. Point W, di W , 

3) The distance from the point W on the diffraction edge to the observer location on the 
ground, Point O, dwo- 
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From these three distances, the difference in source-receiver path length between the direct 
and diffracted sound fields may be computed: 

A = (d| w + dwo) * d[o (7) 

where A > 0 when Point / lies on the wing surface, A = 0 when Point 7 lies on the wing boundary 
edge, and A < 0 when Point I is beyond the wing surface. 


POINT 1 POINT/ 

(Engine Inlet ) (Wing Surface Intersection Point) 



Figure 3. The Point That Is Nearest To The Ray Intersection Point With The Wing Is 
Determined. 



Figure 4. The Differences In Path Length Between The Direct And Diffracted Sound Rays 
Are Used To Calculate The Wing Shielding. 
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From this difference in distances, the Fresnel number is calculated as follows: 


N = 2fi A/Cco 


( 8 ) 


where fi represents the frequency for each 1/3 octave band, in Hz, and Ca> represents the ambient 
speed of sound. 

The attenuation is computed for each 1/3 octave band frequency as follows: 


4 /,)= 


■72 nN 

20 log f— + 5.0 ; N > 0 

tanhT^V 


72^1 N 

20 log - 7 — h== f + 5.0 


tar iy]2n\N 


0 . 


-0.2 < N <0 
N <-0.2 


(9) 


This attenuation is the noise reduction due to a semi-infinite barrier. In this model, 
diffraction around three diffraction edges (wing leading edge, trailing edge, and tip) is included. In 
order to obtain an equivalent total attenuation from the combined effects of the three diffraction 
edges, the individual attenuations for each edge at any frequency fi are combined as follows: 

A roT = -101ogXlCrK /10 ) (10) 


where k = LE, TE, and TIP. 


Wing Reflection Model 

The reflecting panels that represent the wing/flap system are modeled as planar surfaces. The 
geometry definition for each panel consists of the specification of the Cartesian coordinates (x, y, 
z) of each comer point of the panel. The coordinates are defined for the Root Leading Edge (RLE), 
Root Trailing Edge (RTE), Tip Leading Edge (TLE), and Tip Trailing Edge (TTE), as shown in 
Figure 2. The coordinates of the Tip Trailing Edge (TTE) are used only to define the wing 
boundary. Three points define the plane of the wing: RLE, RTE, and TLE. Any twist in the wing 
is ignored. 

The Cartesian coordinates of the wing boundary points are defined relative to the local origin 
of the engine centerline at the exit plane of the fan bypass duct. Local axes are shown in Figure 2. 
Then, the engine fan duct exit and wing coordinates are transformed into a global coordinate 
system consistent with the receiver location on the ground (Point R). This transformation must 
take into account the aircraft attitude and position at the particular time of the observation. 
Therefore, the global coordinates must be calculated at every flyover position. 
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Three wing/flap panels are currently permitted in the model. This allows the description of a 
main wing panel, a flap panel, and a flap tab panel. However, if flaps are not included in the 
model, or are considered to have negligible impact on reflection compared to the wing itself, the 
wing may be modeled using multiple panels, instead. 

The reflection model uses an image source point, which is positioned on a normal to the 
reflecting plane, at a distance equivalent to that between the source and the reflecting plane, but 
located on the opposite side of the reflecting plane. The orientation of the image source is shown 
in Figure 5. 

To obtain the image source point, a normal to the reflecting plane is created, which intersects 
the original source point (S). Then Point IS, the image source, is located on the normal line. The 
points S and IS are equidistant from the reflecting panel. Note that the normal line between Points 
S and IS does not have to intersect the wing panel itself, only the plane extending through the wing 
panel. 



s 


Figure 5. The Orientation of the Image Source Point is Shown Relative to the Original 
Source and the Reflecting Panel. 

First, the components of a vector normal to the reflecting surface are obtained by computing 
the cross-product of vectors along the wing root and leading edge: 


where: 



= 4v ro „, * Az t£ - 

A y LE * 

^ mut 

4 k . 

= *Az le 

-te LE 

*A Z ,) 

root / 


II 

£ 

* 

1 

Ay * 

^ IE 

4v r „„, 


A* mot ~ X RTE 

~ X RLE 



Ay mot = y RTE ' 

~ y RLE 
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root — Zrte z rle 
Ax l£ = X TL£ — X £££ 

= !V TLE ~ y RLE 
AZl£ = Z TLE — Zfiu; 

Next, the coordinates of a point, N, on a line normal to the reflecting surface and passing 
through the engine source point, S, are obtained as follows: 

x N =x s + Ax n 

y N =y s + Ay N 

z N =z s + A z N 

Then, the coordinates of the intersection of the line S-N with the plane of the reflecting 
surface are computed. The intersection point is identified as Point I. The coordinates of the 
intersection point are determined by solving a set of three equations in three unknowns ( x, , y , , 
and z , ). Two equations are produced by the 2-point form of the equation for the line S-N: 


*/-*5 yrzh -Q 
x n ~ x s y n ~ y s 



The other equation comes from the 3-point form of the equation for the reflecting plane: 


X 1 ~~ X RLE 
X RTE ~ X RLE 
X TLE X RLE 


?/ - y RLE 

£/ — Z RLE 

y RTE ~ y RLE 

Zrte ~ Zrle 

yTLE — y RLE 

Ztle ~ Zrle 


Once the coordinates of the intersection point have been determined, the components of the 
vector between the source point, S, and the intersection point, I, may be computed: 


Ax, = x, - x s 

Ay, =y,-y s 
A z, =z,-z s 


Finally, the image source point, IS, may be positioned on the normal to the reflecting surface, 
through the source point, S, by adding the above vector components to the coordinates of the Point 
I: 


x rs = x, + Ax, 
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y ls =yi+ A y, 

z, s = z, +Az ; 

Once the image source is located, it is necessary to determine if a reflection will actually 
occur for the reflecting panel under consideration. A reflection will occur if the line passing 
between the receiver and the image source at a given flyover location intersects the plane of the 
reflecting panel within the boundaries of the actual panel, as shown in Figure 6. The orientation of 
image source and receiver changes at each interval during the flyover, and therefore must be 
recomputed for each flyover position. 



Figure 6. To Determine if a Reflection Will Occur, the Location of the Reflection Point W 
Must Be Determined Relative to the Wing Panel. 

To determine if a reflection will occur, the coordinates of the reflection point (W) on the 
plane of the wing panel must first be obtained, by determining the intersection point of the line 
formed by 1S-R and the plane of the wing panel. The coordinates of the intersection point, W, are 
determined by solving a set of three equations in three unknowns {x w ,y w , and z w ). Two of the 
equations are produced by the 2-point form of the equation for the line IS-R: 


*w- x r >v ~ y* _ 0 

x is ~ x k y is ~ y r 



The other equation comes from the 3-point form of the equation for the reflecting plane: 
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X W X RLE 


yw ~ y RLE 
y RTE ~ y RLE 
yTLE ~ y RLE 


z w z rle 


X RTE X RLE 
X TLE ~ X RLE 


Z RTE ‘“RLE 
Z TLE ~ Z RLE 


= 0 


After the coordinates of the reflection point, W, have been determined, the point must be 
tested to establish whether or not it lies within the boundaries of the reflecting plane. To 
accomplish this, the point on each wing boundary that is nearest to Point W must be located. Each 
of these points (Points Wle, Wte, and W Tn >) is computed by solving a set of three equations in three 
unknowns (e.g., x W££ , y W££ , and z W l£ for the leading edge boundary). The equations are 

obtained by imposing the following conditions, e.g. for the leading edge boundary: 

The line W-Wle must be perpendicular to the wing boundary. This condition is represented 
by setting the dot product of the line W-Wle vector and the wing boundary line vector equal to 
zero, e.g.: 


(c'-VV X W le \ X RLE X TLE )"*" Ov y^iE \y RLE X TLE ) Z W le \ Z RLE Z TLE ) ® 

The point Wle must lie on the wing boundary. This condition is met when the coordinates of 
the point Wle satisfy the 2-point equation of the line representing the wing boundary edge, e.g.: 

x^LE ~ X RLE ywur ~ yRLE 

X TLE ~ X RLE yTLE ~ yRLE 

X W LE ~ X RLE Z W LE -ZRLE 

X TLE ~ X RLE Z TLE ~ Z RLE 

The coordinates of the reflection point, W, are then tested against the coordinates of the three 
nearest boundary points, to determine if Point W is actually located on the wing or flap reflecting 
surface. 

If the reflection point W is contained within the wing panel boundaries, the reflection 
computation may continue. If Point W is outside the wing panel boundaries, no contribution to 
reflected noise occurs, and the algorithm proceeds to consider the next reflecting panel. 

If a reflection has occurred for a particular wing panel, the path length of the reflected ray 
must be determined. The path length is the combined distance from the source to the reflection 
point, S-W, and from the reflection point to the receiver, W-R, as shown in Figure 7. This is 
equivalent to the distance from the image source point, IS, to the receiver point, R. 
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Figure 7. The Path Length of the Reflected Ray is Equivalent to the Distance Between the 
Receiver, R, and the Image Source, IS. 

Thus, the reflected ray path length may be computed from: 

d r =/{ Ax r ) 2 + (Ay,. ) 2 + (Az, ) 2 

where: 

AX r — X IS ~ X R 

Ay, = y, s - y R 

= Z IS ~ Z R 

The strength of the point source, S, is a function of the directivity angle of the ray radiating 
from it. As illustrated in Figure 8, the directivity angles of the direct ray and the reflected ray could 
be different. Therefore, the source strength of the reflected ray must be computed separately from 
that of the direct ray. 



Figure 8. The Source Strength of the Reflected Ray is a Function of the Directivity Angle, 
Which Can Differ From That of the Direct Ray. 


1-12 



To obtain the source strength, the directivity angle of the reflected ray must be computed in 
three-dimensional space. The directivity angle is equivalent to the 3D included angle between the 
engine exhaust centerline and the reflected ray. Then, the reflected ray source strength is 
interpolated from an array of source strength vs. directivity angle. 

To obtain the strength at the ground receiver, the source strength of the reflected ray must be 
corrected for the same propagation effects as the directly-radiated ray, i.e. atmospheric attenuation, 
spherical spreading, and ground reflection. These are all functions of the path length of the 
reflected ray, and thus will have different values than for the direct ray. 

After the corrected strength of each of the reflected rays has been computed, they must be 
combined with the strength of the direct ray, in order to obtain the total Sound Pressure Level 
(SPL) at the receiver. 

The combined contribution of rays from each of the reflecting panels is computed using the 
following equation, for every 1 /3-octave frequency band, k: 

NPANEL . , 

Sf’k.W.TOT = lOIOg X 'O^ 1 * 

I PANELS 

The total contribution from reflected rays may then be combined with the direct ray 
contribution, to yield the total SPL at the receiver. For every 1/3-octave band, the total SPL is 
given by: 

SPL kTOT = 101og[l0 (5ww ' o) + io (s ' Wror/,o) ] 

Noise Prediction 


The method for the preparation of the output noise results for the wing shielding model is as 
follows: 

1. Obtain the geometry and received spectra data from the input files. 

2. For each reception time value, calculate the point / and determine if the ray intersects 
the wing. 

3. Compute the shielding attenuation using equation (10) at the desired values of frequency. 

4. Apply the attenuation to the appropriate value of the received mean-square pressure. 

The output values are the attenuated mean-square acoustic pressure values as a function of 
frequency, reception time, and observer position. 

The method for the preparation of the output noise results for the wing reflection model is as 
follows: 

1. Define the geometry for each reflecting panel (wing, flap, flap tab) 

2. Create the image source for each reflecting panel, at each flyover location 

3. Determine if a reflection will occur for a given panel 
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4. Compute the path length of the reflected ray, for a given panel 

5. Obtain the source strength for the reflected ray, based on directivity 

6. Compute the strength at the receiver, by applying appropriate corrections 

7. Combine the contributions from the reflected ray(s) and the direct ray 

REFERENCES 

1. Beranek, L.L., Noise and Vibration Control, McGraw-Hill Book Company, 1971, pp. 174- 

180. 

2. Maekawa, Z., "Noise Reduction By Screens," Applied Acoustics, Elsevier Publishing Co., 
Ltd., 1968, pp. 157-173. 
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APPENDIX II 


ANOPP 

WING GEOMETRIC EFFECTS MODULE 
USER’S MANUAL 


Il-i 




PURPOSE - WING TAKES NOISE DATA WHICH HAS BEEN PROPAGATED TO THE 
OBSERVER BY THE PROPAGATION (PRO) MODULE 
AND APPLIES CORRECTIONS FOR WING GEOMETRY EFFECTS. 

AUTHOR - DSW(L03/ 02/11) 


INPUT 

USER PARAMETERS 

METHOD OPTION FLAG FOR METHOD TO BE APPLIED 

(NOTE: ORIGIN OF LOCAL COORDINATE SYSTEM FOR 

WING BOUNDARY POINTS IS DIFFERENT FOR 
WING SHIELDING AND WING REFLECTION. 

REFER TO THEORETICAL MANUAL.) 

=1 WING SHIELDING AND DIFFRACTION 
=2 WING REFLECTION 
I PRINT OUTPUT PRINT OPTION CODE 

=0 NO PRINTED OUTPUT 
=1 PRINT INPUT DATA ONLY 
=2 PRINT OUTPUT DATA ONLY 

=3 PRINT BOTH INPUT AND OUTPUT DATA (DEFAULT) 

I UNITS =2HSI , INPUTS ARE IN SI UNITS (DEFAULT) 

=7HENGLISH, INPUTS ARE IN ENGLISH UNITS 
ROOTLE (3) THREE ELEMENT PARAMETER WITH THE X,Y,Z 

COORDINATES OF THE WING ROOT LEADING EDGE (3RS) 
ROOTTE ( 3 ) THREE ELEMENT PARAMETER WITH THE X, Y, Z 

COORDINATES OF THE WING ROOT TRAILING EDGE (3RS) 
TIPLE ( 3 ) THREE ELEMENT PARAMETER WITH THE X,Y,Z 

COORDINATES OF THE WING TIP LEADING EDGE ( 3RS) 

T I PTE (3) THREE ELEMENT PARAMETER WITH THE X,Y,Z 

COORDINATES OF THE WING TIP TRAILING EDGE (3RS) 


THE FOLLOWING PARAMETERS ARE USED FOR METHOD=2 ONLY 


NPANEL 
FLPRLE ( 3 ) 
FLPRTE ( 3 ) 
FLPTLE ( 3 ) 
FLPTTE ( 3 ) 
TABRLE ( 3 ) 
TABRTE ( 3 ) 
TABTLE ( 3 ) 
TABTTE ( 3 ) 


NUMBER OF WING PANELS USED IN WING REFLECTION 


MODEL (I) 

THREE ELEMENT PARAMETER 
COORDINATES OF THE FLAP 
THREE ELEMENT PARAMETER 
COORDINATES OF THE FLAP 
THREE ELEMENT PARAMETER 
COORDINATES OF THE FLAP 
THREE ELEMENT PARAMETER 
COORDINATES OF THE FLAP 
THREE ELEMENT PARAMETER 
COORDINATES OF THE FLAP 
THREE ELEMENT PARAMETER 
COORDINATES OF THE FLAP 
THREE ELEMENT PARAMETER 
COORDINATES OF THE FLAP 
THREE ELEMENT PARAMETER 
COORDINATES OF THE FLAP 


WITH THE X, Y, Z 

ROOT LEADING EDGE (3RS) 

WITH THE X, Y, Z 

ROOT TRAILING EDGE ( 3RS) 

WITH THE X, Y, Z 

TIP LEADING EDGE (3RS) 

WITH THE X, Y, Z 

TIP TRAILING EDGE ( 3RS) 

WITH THE X, Y, Z 

TAB ROOT LEADING EDGE (3RS) 

WITH THE X, Y, Z 

TAB ROOT TRAILING EDGE (3RS) 

WITH THE X , Y , Z 

TAB TIP LEADING EDGE (3RS) 

WITH THE X , Y , Z 

TAB TIP TRAILING EDGE ( 3RS ) 



★ 


DATA BASE UNITS AND MEMBERS 

GEO (BODY) GEOMETRY DATA FOR ALL OBSERVERS RELATIVE 

TO THE AIRCRAFT BODY COORDINATE SYSTEM 
SEE DESCRIPTION IN DATA BASE STRUCTURES. 
(SEE MODULE GEO) 

PRO (PRES) DIMENSIONLESS MEAN SQUARE PRESSURE AT THE 

OBSERVER AS A FUNCTION OF FREQUENCY AND 
TIME. (SEE DESCRIPTION IN DATA BASE 
STRUCTURES . ) 


OUTPUT 

USER PARAMETERS 

NERR = . TRUE . , ERROR ENCOUNTERED, PRO 

TERMINATED ABNORMALLY 
=. FALSE., NO ERRORS ENCOUNTERED, PRO 
TERMINATED SUCCESSFULLY 


DATA BASE UNITS AND MEMBERS 

WING (PRES) DIMENSIONLESS MEAN SQUARE PRESSURE AT THE 

OBSERVER, CORRECTED FOR WING EFFECTS 
AS A FUNCTION OF FREQUENCY AND 
TIME. (SEE DESCRIPTION IN DATA BASE 
STRUCTURES . ) 

SCRATCH (XXXNNN) UNIT SCRATCH CONTAINS THE RESULT OF SUMMING 
NOISE TABLES 

DATA BASE STRUCTURES 

THE FORMAT OF GEO (BODY) IS AS FOLLOWS: 


★ 


if 


RECORD WORD 
1 

1 

2 

3 

4 

5 

6 


DESCRIPTION 

RECORD FORMAT IS I,3RS,I,RS 
OBSERVER INDEX FOR FIRST OBSERVER 
X COORDINATE OF OBSERVER 
Y COORDINATE OF OBSERVER 
Z COORDINATE OF OBSERVER 

NUMBER OF RECEPTION TIMES ASSOCIATED WITH 
THIS OBSERVER (ASSUME VALUE IS N) 
OBSERVER'S HEIGHT 


2 RECORD FORMAT IS *RS 

1 

RECEPTION TIMES FOR CURRENT OBSERVER 
INDEX 
N 


RECORDS 3 THROUGH N+2 CONTAIN GEOMETRY DATA FOR EACH 
RECEPTION TIME. RECORD 3 CONTAINS GEOMETRY DATA FOR 
THE FIRST RECEPTION TIME , RECORD 4 FOR THE SECOND 
RECEPTION TIME , . . . RECORD N+2 FOR THE N TH RECEPTION 
TIME. 

RECORD FORMAT IS *RS 
DISTANCE OF SOURCE FROM OBSERVER 
EMISSION TIME , SEC 
DIRECTIVITY ANGLE, DEG 
ELEVATION ANGLE, DEG 
AZIMUTH ANGLE, DEG 
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4 REPEAT OF RECORD 3 FOR SECOND RECEPTION TIME 


N+3 


RECORD FORMAT IS I,3RS,I,RS 

1 OBSERVER INDEX FOR SECOND OBSERVER 

2 X COORDINATE OF OBSERVER 

3 Y COORDINATE OF OBSERVER 

4 Z COORDINATE OF OBSERVER 

5 NUMBER OF RECEPTION TIMES ASSOCIATED WITH 

THIS OBSERVER (ASSUME VALUE IS M) 


N+4 RECORD FORMAT IS *RS 

1 

RECEPTION TIMES FOR CURRENT OBSERVER 
INDEX 
M 

RECORD N+5 THROUGH RECORD N+M+4 CONTAIN GEOMETRY DATA 
FOR EACH RECEPTION TIME STORED IN THE SAME MANNER AS 
DESCRIBED ABOVE IN RECORDS 3 THROUGH N+2 . 


THE PATTERN AS SEEN IN RECORDS 1 THROUGH N+2 AND RECORDS 
N+3 THROUGH N+M+4 CONTINUES FOR ALL OBSERVERS 

THE FORMAT OF PRO (PRES) AND WING (PRES) IS AS FOLLOWS: 


RECORD 

1 


WORD 

1 


2- (NS+1 ) 


DESCRIPTION 

RECORD FORMAT IS I,*A8 

NUMBER OF NOISE SOURCES PROPAGATED TO 
THE OBSERVERS, NS. 

MODULE NAMES OF NOISE SOURCES PROPAGATED 
TO THE OBSERVERS 


2 


3 


RECORD FORMAT IS 2I,2RS 

1 OBSERVER INDEX FOR THE FIRST OBSERVER 

2 NUMBER OF RECEPTION TIMES ASSOCIATED WITH 
THIS OBSERVER (ASSUME VALUE IS N) 

3 AIR DENSITY AT THE OBSERVER (RE RHO ) 

R 

4 SPEED OF SOUND AT THE OBSERVER (RE C ) 

R 

RECORD FORMAT IS *RS 

1 

RECEPTION TIMES FOR CURRENT OBSERVER 
INDEX 
N 


4 RECORD FORMAT IS *RS 

1 DIMENSIONLESS MEAN SQUARE PRESSURE FOR 
THE FIRST FREQUENCY AND THE FIRST 
RECEPTION TIME 



2 DIMENSIONLESS MEAN SQUARE PRESSURE FOR 
THE SECOND FREQUENCY AND THE FIRST 
RECEPTION TIME 


NF DIMENSIONLESS MEAN SQUARE PRESSURE FOR 
THE LAST FREQUENCY AND THE FIRST 
RECEPTION TIME 


5 


6 


RECORD FORMAT IS *RS 

1 DIMENSIONLESS MEAN SQUARE PRESSURE FOR 

ALL FREQUENCIES FOR THE SECOND 
RECEPTION TIME 
NF 

RECORD FORMAT IS *RS 

1 DIMENSIONLESS MEAN SQUARE PRESSURE FOR 

ALL FREQUENCIES FOR THE THIRD 
RECEPTION TIME 
NF 


N+4 SAME AS RECORD 2 BUT DATA IS FOR SECOND 

OBSERVER 

N+5 SAME AS RECORD 3 BUT DATA IS FOR SECOND 

OBSERVER 

RECORDS 2 THROUGH N+3 REPEAT FOR ALL OBSERVERS. THE 
VALUE OF N DIFFERS FOR EACH OBSERVER. 


ERRORS 

NON -FATAL 

FUNCTIONAL MODULE ERRORS 

1. REQUIRED UNIT MEMBER NOT AVAILABLE 

2. INSUFFICIENT LDS DYNAMIC STORAGE 

3. UNIT MEMBER NOT OF CORRECT FORMAT 

4 . MEMBER MANAGER ERROR OCCURRED ON READING OR OPENING 
A UNIT MEMBER 

7 . ERROR ENCOUNTERED IN BUILDING A UNIT MEMBER 
FATAL - NONE 
LDS REQUIREMENTS 

LENGTH = 9*NFREQ + (2+NFREQ) *NTIME 

WHERE NFREQ = NUMBER OF FREQUENCIES 

NTIME = LARGEST NUMBER OF OBSERVER TIMES 


GDS REQUIREMENTS - NONE 



APPENDIX III 
ANOPP 

WING GEOMETRIC EFFECTS MODULE 
TEST CASE INPUT AND OUTPUT 
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ANOPP JECHO= .FALSE. JLOG=. FALSE. NLPPM=60 $ 

STARTCS $ 

SETSYS JECHO= . FALSE . $ 

$ 

$ 

$ THIS JOB COMPUTES THE CERTIFICATION NOISE LEVELS FOR THE 1992 
$ AST TECHNOLOGY BASELINE BUSINESS JET. THE INPUT 
$ DECK IS SET UP TO TAKE INPUT PARAMETERS THAT MATCH THE INPUT 
$ TO THE GASP PROGRAM TO MAKE IT EASY TO TRANSFER DATA FROM THE 
$ GASP INPUT DECK TO THE ANOPP JOB STREAM. FURTHER EXPLANATION 
$ OF HOW THIS WORKS WILL BE PROVIDED AS THE DATA ARE ENTERED. 

$ 

$ 

$ THE FIRST STEP IS TO ENTER THE GASP NAMLIST DATA. EVERY EFFORT 
$ IS MADE TO KEEP THE DATA IN CONSISTENT FORMAT WITH GASP. 

$ 

$ 

$ NAMELIST "CONT" IS ENTERED FIRST 

$ 

$ 

PARAM IFAA =1 $ CURRENTLY, ONLY OPTIONS 1-4 (APPROACH, 

$ TAKEOFF, SIDELINE, AND LEVEL FLIGHT) ARE 
$ VALID OPTIONS 

PARAM ISI =0 $ SELECT ENGLISH OR SI UNITS 

$ 

$ 

$ NOW, NAMELIST "ENV" IS ENTERED 

$ 

$ 

PARAM TAMB =536.67 $ AMBIENT TEMPERATURE, DEG R 

PARAM PAMB = 2116.22 $ AMBIENT PRESSURE, PSF 

PARAM RH =70. $ RELATIVE HUMIDITY, PERCENT 

PARAM DIST = 100. $ DISTANCE FOR STATIC PREDICTIONS, FT 

$ 

$ 

$ THE ANGLE ARRAY IS NOT ENTERED AS A USER PARAMETER, BUT AS A UNIT 
$ MEMBER (FILE) ALONG WITH THE DESIRED FREQUENCIES AS FOLLOWS: 

$ 

$ 

UPDATE NEWU=SFIELD SOURCE=* $ 

-ADDR OLDM=* NEWM=FREQ FORMAT=4H*RS$ $ 1/3 OCTAVE CENTER FREQUENCIES 

50. 63. 80. 100. 

125. 160. 200. 250. 315. 400. 500. 630. 800. 1000. 

1250. 1600. 2000. 2500. 3150. 4000. 5000. 6300. 8000. 10000. $ 

-ADDR OLDM= * NEWM=THETA FORMAT=4H*RS$ $ POLAR DIRECTIVITY ANGLES 
10. 20. 30. 40. 50. 60. 70. 80. 90. 100. 110. 120. 130. 140. 

150. 160. $ 

-ADDR OLDM=* NEWM=PHI FORMAT=4H*RS$ $ AZIMUTH DIRECTIVITY ANGLES 
0. $ SOURCES ARE AXISYMMETRIC 

END* $ 

$ 

$ 

$ IN ADDITION, THE TEMPERATURE AND RELATIVE HUMIDITY MUST BE ENTERED 
$ AS A UNIT MEMBER (FILE) BECAUSE ANOPP ASSUMES YOU ALWAYS WANT TO 
$ USE A LAYERED ATMOSPHERE 
$ 
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$ 

UPDATE NEWU=ATM SOURCE=* $ 

-ADDR OLDM=* NEWM=IN FORMAT=4H3RS$ $ 

0. 536.67 70. $ ALTITUDE, TEMPERATURE, RELATIVE HUMIDITY 

END* $ (ONLY ONE RECORD IS NEEDED FOR UNIFORM ATMOSPHERE) 

$ 

$ 

$ THE NAMELIST VARIABLES FOR "SYS" ARE ENTERED NEXT 


$ 

$ 

PARAM 

NTYE 


1 

$ 

ONLY CURRENT OPTION IS TURBOJET (OR MIXED 





$ 

STREAM TURBOFAN) 

PARAM 

ICOMP 

= 

1 , 4,5 

$ 

ONLY CURRENT OPTIONS ARE FAN, COMBUSTOR, 





$ 

OR JET 

PARAM 

ENP 

= 

2 . 

$ 

NUMBER OF ENGINES 

PARAM 

ANENGI 

= 

0 . 

$ 

ANGLE BETWEEN ENGINE INLET AND AIRCRAFT, 





$ 

DEGREES 

PARAM 

ANENGE 

= 

0 . 

$ 

ANGLE BETWEEN ENGINE EXHAUST AND AIRCRAFT 





$ 

DEGREES 

PARAM 

WGMAX 


28700. 

$ 

AIRCRAFT MAX. GROSS WEIGHT AT T/O, LB 

PARAM 

AMACH 

= 

0.2086 

$ 

AIRCRAFT MACH NUMBER (CAN ALSO SPECIFY 

$ 

$ 




$ 

PARAMETER "VEL" IN FT/ SEC) 


$ THE VARIABLES LOCENG, XL, YL, ZL, I PHASE, AND IDOP ARE NOT 
$ APPLICABLE TO THE CURRENT ANOPP MODEL 
$ 

$ 

$ THE NAMELIST "FPRO” IS ENTERED NEXT 


S 

$ 

PARAM 

ID PRO 

= 0 

$ 

$ 

STRAIGHT LINE PROFILE (USER CAN ALSO 
SPECIFY PROFILE USING UNIT MEMBER) 

PARAM 

FLTANG 

= -3.0 

$ 

FLIGHT PATH ANGLE, DEGREES 

PARAM 

ANGAFT 

= 4.2 

$ 

AIRCRAFT ANGLE OF ATTACK, DEGREES 

PARAM 

TOROLL 

= 4921.3 

$ 

LENGTH OF TAKEOFF ROLL, FT 

PARAM 

APDIST 

= 10685. 

$ 

INITIAL AIRCRAFT APPROACH RANGE, FT 

PARAM 

$ 

$ 

XALT 

= 1000. 

$ 

AIRCRAFT ALTITUDE FOR LEVEL FLYOVER, 


$ NOW THE ENGINE THERMODYNAMIC DATA ARE ENTERED. NAMELIST "FAN" FOR 
$ PREDICTING FAN NOISE IS ENTERED FIRST 


$ 

$ 

PARAM 

IGV 


0 

$ 

FAN HAS NO INLET GUIDE VANES 

PARAM 

NBF 

= 

30 

$ 

NUMBER OF FAN BLADES 

PARAM 

NVAN 

= 

61 

$ 

NUMBER OF STATOR VANES 

PARAM 

RSS 

= 

170. 

$ 

ROTOR/ STATOR SPACING IN 

PARAM 

WAFAN 

= 

91.455 

$ 

FAN INLET WEIGHT FLOW, LB/S 

PARAM 

RPM 

= 

6982 . 

$ 

FAN PHYSICAL SPEED, RPM 

PARAM 

FPR 

= 

1.239 

$ 

FAN PRESSURE RATIO 

PARAM 

FANDIA 

- 

2.455 

$ 

FAN DIAMETER, FT 

PARAM 

TIPMD 

= 

1.446 

$ 

FAN TIP MACH NUMBER AT DESIGN POINT 

PARAM 

$ 

$ 

FANEFF 


0.8104 

$ 

FAN EFFICIENCY 


$ NOW NAMELIST "BURNER" FOR THE COMBUSTOR 


m-2 



$ 

$ 

PARAM WACOMB =14.310 $ COMBUSTOR WEIGHT FLOW, LB/S 

PARAM T3 = 1002.9 $ COMBUSTOR INLET TEMPERATURE, DEG R 

PARAM T4 = 2074.9 $ COMBUSTOR EXIT TEMPERATURE, DEG R 

PARAM P3 = 13737.3 $ COMBUSTOR INLET TOTAL PRESSURE, PSF 

$ 

$ 

$ AND FINALLY "JET" FOR THE MIXED STREAM JET NOISE 

$ 

$ 

PARAM VJ =692.4 $ FULLY EXPANDED JET VELOCITY, FPS 

PARAM TJ =754.3 $ JET TOTAL TEMPERATURE, DEG R 

PARAM DJ =1.8507 $ JET OUTER DIAMETER, FT 

PARAM GAMJ =1.333 $ JET RATIO OF SPECIFIC HEATS 

$ 

$ 

$ THERE ARE MANY ADDITIONAL PARAMETERS AND NAMELISTS FOR GASP. HOWEVER, 

$ THEY ARE NOT RELEVANT TO ANOPP CAB ABILITIES OR HAVE NOT BEEN IMPLEMENTED 
$ IN THIS JOB STREAM. THE REMAINING STATEMENTS ARE REQUIRED TO CONVERT 
$ TO ANOPP INPUT AND TO EXECUTE THE ANOPP MODULES. THEY NEED NEVER BE 
$ CHANGED UNLESS THE JOB STREAM CAPABILITY IS MODIFIED. 

$ 

$ 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

$ 

$ CONVERSION OF GASP USER PARAMETERS TO ANOPP INPUT 

$ 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

$ 

$ 

$ FIRST, THE OBSERVER POSITIONS CORRESPONDING TO THE FOUR FAR 36 FLIGHT 
$ PROFILES ARE DEFINED 
$ 

$ 

IF ( IFAA .GT. 1 ) GOTO A1 $ 

$ 

$ 

$ APPROACH OBSERVER COORDINATES 

$ 

$ 

UPDATE NEWU=OBSERV SOURCE=* $ 

-ADDR OLDM= * NEWM=COORD FORMAT=4H3RS$ $ 

-7516. 0. 4. $ 

END* $ 

GOTO A4 $ 

$ 

A1 CONTINUE $ 

IF ( IFAA .GT. 2 ) GOTO A2 $ 

$ 

$ 

$ TAKEOFF OBSERVER COORDINATES 

$ 

$ 

UPDATE NEWU=OBSERV SOURCE=* $ 

-ADDR OLDM=* NEWM=COOED FORMAT=4H3RS$ $ 
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21325 . 


0. 4. $ 


END* $ 

GOTO A 4 $ 

$ 

A2 CONTINUE $ 

IF ( IFAA .GT. 3 ) GOTO A3 $ 

$ 

$ 

$ SIDELINE OBSERVER COORDINATES 

$ 

$ 

UPDATE NEWU=OBSERV SOURCE=* $ 


-ADDR OLDM= 

* NEWM= 

: COORE 

> FORMAT=4H3RS$ 


6000. 

1476. 

4. 

$ 


7000. 

1476. 

4. 

$ 


8000. 

1476. 

4. 

$ 


9000. 

1476. 

4. 

$ 


10000. 

1476. 

4. 

$ 


11000. 

1476 . 

4. 

$ 

END* 

$ 




GOTO 

A4 $ 




$ 





A3 CONTINUE 

$ 




$ 

$ 

$ LEVEL FLYOVER OBSERVER COORDINATES 

$ 

$ 

UPDATE NEWU=OBSERV SOURCE=* $ 

-ADDR OLDM= * NEWM=COORD FORMAT=4H3RS$ $ 

0. 0. 4. $ 

END* $ 

A4 CONTINUE $ 

$ 

$ 

$ NOW SOME STANDARD CONTROL PARAMETERS ARE DEFINED 

$ 

$ 


PARAM 

PIE 

= 

3.14159 

$ 

VALUE OF PI 


PARAM 

AE 

= 

1. 


$ 

SET REFERENCE AREA TO ONE SQUARE 

FOOT 

PARAM 

RS 

= 

DIST 


$ 

SET SOURCE RADIUS DISTANCE INPUT 

VALUE 

PARAM 

TA 

= 

TAMB 


$ 

DEFINE AMBIENT TEMPERATURE 


EVALUATE 

RHOA 

= 

PAMB 

/ 

TAMB 

/ 1716.22 







$ 

COMPUTE AMBIENT DENSITY 


EVALUATE 

CA 

= 

1116 . 

22 

* SQRT ( TAMB / 518.67 ) 







$ 

COMPUTE AMBIENT SPEED OF SOUND 


EVALUATE 

NENG 

= 

INT ( 

ENP ) 








$ 

MAKE ENGINE NUMBER INTEGER 


PARAM 

MA 

= 

AMACH 

$ 

DEFINE MACH NUMBER 


PARAM 

IOUT 

= 

1 


$ 

PRINT DB VALUES ONLY 


IF { ISI 

.NE. 0 

) 

GOTO 

B1 

$ 

SET UNITS FLAG 


PARAM 

IUNITS 

= 

7HENGLISH $ 



GOTO B2 

$ 







B1 CONTINUE $ 







PARAM 

IUNITS 

= 

2HSI 

$ 




B2 CONTINUE $ 








$ 
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$ THE ENGINE PARAMETERS ARE CONVERTED TO ANOPP FORM 
$ 

$ FIRST, THE FAN 
$ 

$ 


EVALUATE 

AFAN 


PIE * 

FANDIA* *2 / 4. 

$ COMPUTE FAN REFERENCE AREA 

PARAM 

DIAM 


FANDIA 

$ 

DEFINE FAN DIAMETER 

PARAM 

MD 

= 

TIPMD 


$ 

DEFINE TIP MACH NUMBER AT DESIGN POINT 

EVALUATE 

RSS 

= 

RSS / 

100 . 
$ 

CONVERT ROTOR/STATOR SPACING TO RATIO 

EVALUATE 

MDOT 

= 

WAFAN 

/ 

32 . 
$ 

17 / RHOA / CA 
NORMALIZE WEIGHT FLOW 

EVALUATE 

DELTAT 


( FPR' 

* * 0 

.2857 - 1. ) / FANEFF 

$ COMPUTE FAN TEMPERATURE RISE 

PARAM 

NB 

= 

NBF 


$ 

SET NUMBER OF BLADES 

PARAM 

NV 

= 

NVAN 


$ 

SET NUMBER OF VANES 

EVALUATE 

IGV 

= 

IGV + 

1 

$ 

SET IGV FLAG 

EVALUATE 

N 

= 

{ RPM 

/ 

60. 

$ 

. ) * DIAM / CA 

COMPUTE NORMALIZED ROTATIONAL SPEED 

PARAM 

INCT 

= 

. FALSE . 

$ 

TURN OFF COMBINATION TONES 


$ 

$ 

$ NOW THE COMBUSTOR 

$ 

$ 

EVALUATE A =0.1* AFAN 

$ ARBITRARY AREA DEFINED 
EVALUATE MDOTC = WACOMB / 32.17 / RHOA / CA 

$ WEIGHT FLOW NORMALIZED (NOTE: ANOPP USES 
$ “ MDOT " FOR BOTH THE FAN AND COMBUSTOR - 

$ COMBUSTOR MASS FLOW IS RENAMED TO AVOID 
$ OVERWRITE 

EVALUATE PI = P3 / PAMB 

$ NORMALIZE INPUT PRESSURE 

EVALUATE TI = T3 / TAMB 

$ NORMALIZE INPUT TEMPERATURE 

EVALUATE TCJ = T4 / TAMB 

$ NORMALIZE OUTPUT TEMPERATURE 
PARAM TDDELT =1.0 $ USE THIS PARAMETER - SET TO 1 

$ 

$ 

$ JET PARAMETERS 

$ 

$ 

EVALUATE AJ = PIE * DJ ** 2 / 4 . 

$ COMPUTE JET AREA 

EVALUATE TJ = TJ / TAMB 

$ NORMALIZE JET TOTAL TEMPERATURE 
EVALUATE VJ = VJ / CA $ NORMALIZE JET VELOCITY 

EVALUATE RHOJ =1. / ( TJ - { GAMJ - 1 ) / 2. * VJ**2 ) 

$ COMPUTE NORMALIZED JET DENSITY 

PARAM CIRCLE = .TRUE. $ REQUEST SINGLE JET FORM STONE'S METHOD 

$ 

$ 

$ LOAD UNITS FROM DATA LIBRARY 
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$ 

$ 

LOAD /LIBRARY/ SAE PROCLIB STNTBL $ 

$ 

$ 

$ PREDICT SOURCE NOISE 

$ 

$ 

PARAM IDBB = .FALSE. $ 

PARAM IDRS = .FALSE. $ 

EXECUTE HDNFAN HDNFAN=FANIN $ 

PARAM IDBB = .TRUE. $ 

PARAM IDRS = .TRUE. $ 

PARAM INRS = .FALSE. $ 

PARAM INDIS = .FALSE. $ 

PARAM INBB = .FALSE. $ 

EXECUTE HDNFAN HDNFAN= FANOUT $ 

EVALUATE RS = DIST * SQRT ( 10. ) $ 

EXECUTE GECOR MDOT=MDOTC $ 

PARAM RS = DIST $ 

EXECUTE SGLJET $ 

$$$$$ EXECUTE STNJET A1=AJ DE1=DJ DH1=DJ V1=VJ T1=TJ RH01=RH0J $ 

$ 

$ 

$ NOW, THE ATMOSPHERIC CONDITIONS AND THE ATMOSPHERIC ABSORPTION 
$ COEFFICIENTS ARE COMPUTED 
$ 

$ 

EXECUTE ATM P1=PAMB $ 

EXECUTE ABS $ 

$ 

$ 

$ THE FLIGHT PATH AND GEOMETRY IS NOW DEFINED 

$ 

$ 

EVALUATE VA = AMACH * CA $ DEFINE AIRCRAFT SPEED 

IF ( IFAA .GT. 1 ) GOTO Cl $ 

EVALUATE XA = 0. - APDIST 

$ DEFINE STARTING DISTANCE FOR APPROACH 
EVALUATE ZA = - APDIST * SIN ( FLTANG ) 

$ DEFINE ALTITUDE AT BEGINNING OF APPROACH 

EVALUATE THW = 0. - FLTANG 

$ DEFINE FLIGHT PATH ANGLE 
PARAM PLG = 4HDOWN $ LANDING GEAR IS DOWN 

PARAM TLG = -1. $ LANDING GEAR CHANGED BEFORE START 

PARAM JF =200 $ ALLOW 200 TIME STEPS 

PARAM ZF =0. $ STOP AT TOUCHDOWN 

PARAM START =0. $ START EPNL CALCULATION 

PARAM STOP =80. $ STOP EPNL CALCULATION 

GOTO C3 $ 

Cl CONTINUE $ 

IF ( IFAA .GT. 3 ) GOTO C2 $ 

PARAM XA = TOROLL 

$ DEFINE STARTING DISTANCE FOR TAKEOFF 


PARAM ZA = 0. 


DEFINE ALTITUDE AT BEGINNING OF TAKEOFF 


PARAM THW = FLTANG 
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$ DEFINE FLIGHT PATH ANGLE 


PARAM 

PLG 

= 

4HUP 


$ LANDING GEAR IS UP 

PARAM 

TLG 


- 1 . 


$ LANDING GEAR CHANGED BEFORE START 

PARAM 

JF 

= 

200 


$ ALLOW 200 TIME STEPS 

PARAM 

ZF 

- 

32000. 


$ STOP AT 32000. FT 

PARAM 

XF 

= 

32000. 


$ 

GOTO C3 : 

$ 





C2 CONTINUE $ 





EVALUATE 

XA 

= 

XALT * SIN 

( 

-5. ) 





$ 

DEFINE STARTING DISTANCE FOR LEVEL FLIGHT 

PARAM 

ZA 

= 

XALT 

$ 

DEFINE ALTITUDE FOR LEVEL FLYOVER 

PARAM 

THW 

= 

0 . 

$ 

DEFINE FLIGHT PATH ANGLE 

PARAM 

PLG 

= 

4HUP 

$ 

LANDING GEAR IS UP 

PARAM 

TLG 

= 

-100. 

$ 

LANDING GEAR CHANGED BEFORE START 

PARAM 

JF 


200 

$ 

ALLOW 200 TIME STEPS 

EVALUATE 

XF 

= 

0. - XA 

$ 

STOP AT SAME DISTANCE FROM MIC. 

PARAM 

START 


-9999. 

$ 

SET ARBITRARY START TIME 

C3 CONTINUE $ 





PARAM 

ALPHA 

= 

ANGAFT 

$ 

SET ANGLE OF ATTACK 

PARAM 

ENGNAM 

= 

3HXXX 

$ 

SET DEFAULT ENGINE NAME IN SFO 

PARAM 

ICOORD 

= 

1 

$ 

REQUEST BODY AXIS 


$ 

$ 

$ NOW GENERATE GEOMETRY 

$ 

$ 

PARAM I PRINT = 1 $ PRINT INPUT ONLY 

$ 

EXECUTE SFO VI=VA XI=XA ZI=ZA VF=VA $ 
EXECUTE GEO $ 

$ 

$ 

$ NOW ENTER PROPAGATION PARAMETERS 

$ 

$ 


PARAM 

NCOMP 

= 

1 

$ 

NUMBER OF NOISE COMPONENTS TO BE PROPAGATED 

PARAM 

ABSORP 

= 

. TRUE . 

$ 

INCLUDE ABSORPTION 

PARAM 

GROUND 

= 

. TRUE . 

$ 

INCLUDE GROUND EFFECTS 

PARAM 

PROSUM 

= 

6HFANOUT 






$ 

FOUR NOISE SOURCES 

PARAM 

IOSPL 

= 

. TRUE . 

$ 

COMPUTE OVERALL SPL 

PARAM 

I AWT 

= 

. TRUE . 

$ 

COMPUTE A-WEIGHTED OASPL 

PARAM 

IPNL 

= 

. TRUE . 

$ 

COMPUTE PNL 

PARAM 

PROPRT 

= 

1 

$ 

ONLY PRINT PROPAGATION MODULE INPUT 

PARAM 

LEVPRT 

= 

1 

$ 

ONLY PRINT NOISE LEVELS MODULE INPUT 

PARAM 

EFFPRT 

= 

1 

$ 

ONLY PRINT EFFECTIVE NOISE MODULE INPUT 


$ 

$ 

$ UNIT FLI MUST BE MODIFIED TO SET ONLY ONE SOURCE TIME 

$ 

$ 

UPDATE NEWU=FLIMOD OLDU=FLI ALL SOURCE=* $ 

-OMIT FLIXXX $ 

-ADDR OLDM=* NEWM=FLIXXX FORMAT=l 1H6RS , A4 , 2RS$ $ 

0. 0.2 1. 1116. .00238 .1 4 HUP 0. 0. $ 


HI-7 



END* $ 

$ 

$ 

$ AND CALL PROPAGATION MODULE 
$ 

$ 

$ 

EXECUTE PRO GEOM=BODY FLI=FLIMOD $ 
$ 


PARAM 

ROOTLE 

= 

4.6, 

GO 

T* 

1 

-17.1 

$ 

WING 

ROOT LEADING EDGE 

PARAM 

ROOTTE 

= 

-8.2, 

-4.3, 

-17.1 

$ 

WING 

ROOT TRAILING EDGE 

PARAM 

TIPLE 

= 

-7.4, 

-4.3, 

21.0 

$ 

WING 

TIP 

LEADING EDGE 

PARAM 

TIPTE 

= 

-12.6, 

-4.3, 

21.0 

$ 

WING 

TIP 

TRAILING EDGE 

PARAM 

FLPRLE 

= 

-7.7, 

-3.9, 

-16.9 

$ 

FLAP 

ROOT LEADING EDGE 

PARAM 

FLPRTE 

= 

-10.1, 

-2.5, 

-16.9 

$ 

FLAP 

ROOT TRAILING EDGE 

PARAM 

FLPTLE 

= 

-11.2, 

-4.0, 

11.4 

$ 

FLAP 

TIP 

LEADING EDGE 

PARAM 

FLPTTE 

= 

-12 .5, 

-3.3, 

11.6 

$ 

FLAP 

TIP 

TRAILING EDGE 

PARAM 

TABRLE 

= 

-10.3, 

-2.5, 

-16.9 

$ 

FLAP 

TAB 

ROOT LEADING EDGE 

PARAM 

TABRTE 

= 

-11.4, 

-0.9, 

-16.9 

$ 

FLAP 

TAB 

ROOT TRAILING EDGE 

PARAM 

TABTLE 

= 

-12.7, 

-3.2, 

11.6 

$ 

FLAP 

TAB 

TIP LEADING EDGE 

PARAM 

TABTTE 

= 

-13.3, 

-2.4, 

11.8 

$ 

FLAP 

TAB 

TIP TRAILING EDGE 

PARAM 

METHOD 

= 

2 , I PRINT = 

3, NPANEL = 3 

l $ 



$ 

EXECUTE WING $ 

$ 

$ MEMLIST LIST=FULL WING (PRES) FORMAT= 4 H * RS $ $ 
$ PARAM MEMSUM=4HPRO 4HPRES $ 

$ EXECUTE LEV $ 

PARAM MEMSUM=4HWING 4HPRES $ 

EXECUTE LEV $ 

$ 

$ END OF JOB 

$ 

$ 

ENDCS $ 
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TIPLE = -7.40000 -4.30000 21.00000 
TIPTE = -12.60000 -4.30000 21.00000 
FLPRLE = -7.70000 -3.90000 -16.90000 
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OB S. TIME WING GEOMETRY SPL CORRECTION FOR EACH FREQUENCY 

21.20 SPL .00 .00 .00 .00 .00 .00 .00 .00 

SPL .00 .00 .00 .00 .00 .00 .00 .00 

SPL .00 .00 .00 .00 .00 .00 .00 .00 
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OBS. TIME WING GEOMETRY SPL CORRECTION FOR EACH FREQUENCY 

24.70 SPL .00 .00 .00 .00 .00 .00 .00 .00 

SPL .00 .00 .00 .00 .00 .00 .00 .00 

SPL .00 .00 .00 .00 .00 .00 .00 .00 



OB S. TIME WING GEOMETRY SPL CORRECTION FOR EACH FREQUENCY 

25.20 SPL .00 .00 .00 .00 .00 .00 .00 .00 

SPL .00 .00 .00 .00 .00 .00 .00 .00 

SPL .00 .00 .00 .00 .00 .00 .00 .00 
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OBS. TIME WING GEOMETRY SPL CORRECTION FOR EACH FREQUENCY 

33.20 SPL .00 .00 .00 .00 .00 .00 .00 .00 

SPL .00 .00 .00 .00 .00 .00 .00 .00 

SPL .00 .00 .00 .00 .00 .00 .00 .00 



OB S. TIME WING GEOMETRY SPL CORRECTION FOR EACH FREQUENCY 

33.70 SPL .00 .00 .00 .00 .00 .00 .00 .00 

SPL .00 .00 .00 .00 .00 .00 .00 .00 

SPL .00 .00 .00 .00 .00 .00 .00 .00 
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NOISE SOURCES USED IN CALCULATING THE NOISE LEVELS 
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